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The  purpose  of  this  report  is  to:  (1)  document  the  0^(a*A)  diagnostics  and 

calibration  procedures  used  for  the  ROTOCOIL  test  progreun,  (2)  establish  a 

baseline  calibration  apparatus  and  procedure  for  various  COIL  research  groups, 

(3)  review  the  status  of  absolute  0  (a*A)  detection  methods,  and  (4)  introduce 

2 

spectroscopic  and  radiometric  techniques  to  the  novice  researcher. 
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1.0  INTRODUCTION 

The  Chemical  Oxygen-Iodine  Laser  (COIL)  is  the  shortest  wavelength  high  energy 
chemical  laser  in  existence  today  and  has  great  potential  for  Air  Force  direc¬ 
ted  energy  weapons  applications.  The  Phillips  Laboratory  (PL),  Kirtland  Air 
Force  Base,  New  Mexico,  is  vigorously  pursuing  the  development  of  efficient, 
scalable  COIL  devices  with  high  output  power  and  excellent  beam  quality. 

While  chemical  oxygen-iodine  lasers  are  in  the  early  stages  of  large  scale 
development,  a  25-kU  device  is  being  tested  at  the  PL  auid  scaling  to  larger 
devices  is  under  active  investigation. 

The  p>erformance  of  COIL  devices  strongly  depends  on  the  absolute  concentration 

and  yield  of  the  metastable  energy  carrier,  singlet  oxygen.  The  total  energy 

density  available  in  the  COIL  flow  is  determined  by  the  concentration  of 

O^Ia^A);  however,  the  extractable  energy  is  limited  by  the  yield,  defined  as 

the  ratio  [0  (^A)]/[0  (}C^£)]  (Ref.  1).  Since  typical  chemical  oxygen  genera- 
2  2 

tors  operate  near  the  threshold  yield,  the  laser  power  is  sensitive  to  small 
changes  in  singlet  oxygen  concentrations  (Ref.  2).  Without  reliable  measure¬ 
ments  of  singlet  oxygen,  it  is  Impossible  to  separate  the  effects  of  oxygen 
generator  performance,  iodine  n(.2zle  mixing  efficiency,  and  cavity  transport 
on  attained  laser  power.  The  performance  of  chemical  singlet-oxygen  genera¬ 
tors  is  best  determined  by  the  absolute  measurement  of  singlet-oxygen  concen¬ 
trations.  Additionally,  the  study  of  COIL  gas  phase  kinetics  often  requires 

diagnostics  for  absolute  0  (^A)  concentrations  (Ref.  3).  Thus,  the  accurate, 

^  1 

absolute  measurement  of  0^(a  A)  is  critical  to  predicting,  characterizing,  and 
modeling  the  performance  of  COIL  devices. 

One  method  to  measure  O^(a^A)  concentration  is  to  observe  its  spontaneous 
emission  with  a  detector  calibrated  using  Electron  Paramagnetic  Resonance 
(EPR)  spectroscopy  (Ref.  4).  Both  ground  state,  X^E,  atnd  excited  state,  a^A 
oxygen  have  been  studied  using  EPR.  In  fact,  ground  state  oxygen  is  often 
used  as  a  quantitative  standard  for  EPR  (Ref.  5).  The  EPR  spectrum  of  O^(^A) 
was  first  reported  in  1965  (Ref.  6),  and  has  received  considerable  attention 
both  theoretically  and  experimentally  (Refs.  7-9).  Indeed,  EPR  spectroscopy 
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is  currently  the  nost  reliable  technique  for  aeasurlng  the  concentration  of 

O^A)  under  low  pressure,  gas  phase  conditions  (Ref.  4).  By  observing  the 

relative  EPR  spectral  areas  of  O^(X^)  and  O^(a^A)  In  a  pure  oxygen  flow,  a 

standard  for  the  absolute  concentrations  of  these  species  Is  established 

(Ref. 4).  If  the  infrared  emission  from  the  0  (a-X)  transition  Is  slmul- 

2 

taneously  observed,  then  a  method  for  calibrating  the  emission  detector 
Is  provided.  This  technique  has  been  used  extensively  to  characterize  the 
performance  of  CX)IL  devices  (Refs.  4,  10-12). 

Conceptually  the  method  is  straightforward  but  many  experimental  and  analyti 
cal  issues  limit  the  accuracy.  In  this  report,  an  attempt  will  be  made 
to  clearly  define  and  reduce  the  most  Important  systematic  and  statistical 
errors  and  to  provide  suggestions  for  designing  a  reliable  experimental 
apparatus  for  the  calibration  and  operation  of  such  a  diagnostic  on  COIL 
devices.  Specifically,  the  following  issues  will  be  addressed:  derivation 
of  the  theoretical  expressions  required  for  absolute  0  (a*A)  concentration 
measurements  (Section  2.0),  characterization  of  a  reliable,  accurate  callbra 
tion  apparatus  (Section  3.0),  analysis  of  observed  spectra  and  emission 
signals  (Section  4.0),  and  a  detailed  error  analysis  (Section  5.0). 
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2.0  BACKGROUND  THEORY 

Theoretical  expressions  relating  measured  EPR  signals  to  absolute  oxygen  con¬ 
centrations  and  observed  infrared  detector  signals  will  be  reviewed  in  this 
section.  A  background  discussion  of  spectroscopy,  kinetics,  and  radioaetry 
is  provided  to  assist  in  this  development. 

2.1  ELECTRCMIIC  STRUCTURE  OF  OXYGEN 


2  2  4 

Oxygen  atoms  have  ground  state  electronic  configurations  Is  2s  2p  .  The  two 
unfilled  2p  orbitals  are  equivalent  to  two  electrons  with  J  »  1,  s  «  1/2 
yielding  IS  distinct  states  (Ref.  13).  Using  a  Slater  Diagram,  the  three  mul- 
tiplets  V,  ^D,  and  ^S  are  readily  identified.  The  spin-orbit  splitting  of 
the  ^  term  to  form  the  and  levels  is  small  with  the  energies 

of  these  levels  being  0,  158.5  and  221.5  cm"^,  respectively  (Ref.  14)  .  The 
triplet  term,  is  the  lowest  in  energy  with  the  singlet  ^D  and  ^S  terms 
higher  by  1.967  and  4.188  eV,  respectively  (Ref.  14). 

The  molecular  orbitals  for  diatomic  oxygen  may  be  formed  from  linear  combina¬ 
tions  of  the  atomic  orbitals  (Refs.  10-16).  Figure  1  Illustrates  the  atomic 
orbitals  and  the  resulting  molecular  orbitals  for  oxygen.  The  filling  of 

these  molecular  orbitals  by  the  16  electrons  of  0.  is  denoted  by  the  arrows 

2 

indicating  spin  up  (m^  •  *i/2)  or  down  (m^  ■  -1/2).  The  configuration  of 
Figure  1  may  be  written 

(le  )*(lir  )*(2<r  )^(2<r  )^(3«r  )*(1k  )*(la  )*(3<r  )° 
gugugugu 

where  the  superscripts  denote  the  number  of  electrons  in  the  given  orbital. 

Typically  this  notation  is  shortened  to  2420  for  the  occupancy  of  the  last 

four  orbitals.  The  antibonding  orbital,  is  ,  is  only  partially  filled  and 

9 

there  are  six  unique  ways  for  the  two  electrons  to  occupy  this  orbital.  The  6 
arrangements  resulting  in  6  possible  electronic  levels  are  shown  in  Figure  2. 
The  ai^ropriate  antisymmetric  wavefunctions  are  described  in  Reference  17. 

The  ground  electronic  state  is  '  and  the  first  two  excited  states  are  a^A 

9  9 
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Figure  1 


Figure  2 
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Molecular  orbitals  for  0  from  linear  combinations  of  atomic 

2 

orbitals.  Arrows  indicate  occupancy  by  am  electron  w’th  m^  =  .^1/2. 
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Electronic  states  for  molecular  oxygen  and  the  corresponding 
electronic  configuration  within  the  antibonding  lit  orbital. 
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and  b^£  * .  Potential  energy  curves  for  these  and  higher  lying  states  of 
diatonic  oxygen  are  shovm  in  Figure  3  (Ref.  13). 


Since  the  three  lowest  electronic  terms  of  oxygen,  X,  a,  and  b,  all  arise  from 
the  same  configuration,  electric  dipole  transitions  are  not  allowed  between 
these  terms.  The  O^(a-X)  transition  is  doubly  forbidden  (  g  -»  g  and  singlet 
to  triplet  transitions  are  both  forbidden  for  dipole  transitions)  in  the 
electric  dipole  since  there  is  no  change  in  Inversion  symmetry,  but  a  change 
of  one  in  total  spin,  and  a  change  of  two  in  total  angular  momentum.  As  a 
result,  the  radiative  lifetime  of  O^(a^A)  is  a  64  min  (Ref.  17).  The  Einstein 
A  coefficient  for  the  0  (a^A,  v’  =  0)  to  (X^  ”,  v"  »  0)  transition  as  A^  = 
2.58  X  10~*  1/s  was  measured  with  a  15  percent  accuracy  (Ref.  17).  Radiative 
lifetimes  as  long  as  90  min  have  been  ’Reported  (Ref.  12).  Since  the  radiative 
lifetime  is  so  long,  self  absorption  t.riects  are  negligible  for  optical  path 


lengths  <  100  km  (Ref.  10).  The  first  highly  allowed  electronic  transition  in 

oxygen  occurs  in  the  ultraviolet  between  the  ”  and  X^  ”  states  (Ref.  13). 

3L.  -  a  a 

The  state  arises  from  the  molecular  orbital  configuration  2330  where  one 

electron  has  been  promoted  from  the  lir  bonding  orbital  to  the  la  auntibonding 

«  9 

orbital. 


The  O^(a-X)  and  O^(b-X)  trwsitions  were  recognized  early  as  Important  atmo¬ 
spheric  absorption  features  (Refs.  18-20).  Detailed  absorption  experiments 
(Refs.  18-19)  have  been  conducted  and  Franck-Condon  factors  reported  (Ref. 
21).  Typical  spectra  of  the  O^(a-X)  and  (b-X)  transitions  under  moderate 
resolution  are  shown  in  Figure  4  and  spectroscopic  constants  for  the  X,  a, 
and  b  states  are  reported  in  Table  1  (Ref.  22).  Note  that  the  equilibrium 
intemuclear  separation  is  nearly  identical  for  the  X  and  a  states.  As  a 
result,  the  (0,0)  and  (1,1)  trauisitlons  are  highly  favored  and  a  large  frac¬ 
tion  of  the  total  emission  from  a  vibrationally  thermalized  O^(^A)  state 
occurs  at  1.268  iim.  This  lack  of  dispersion  greatly  assists  in  developing 
emission  diagnostics. 
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Figure  3.  Potential  energy  curves  for  nolecular  oxygen.  (Fro«  Ref.  7) 
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Wavelength  (X) 

a)  Spectrum  of  O^(b-X)  showing  the  v’  =0,  v"  =0  band. 


Wavelength  (X) 


b)  Spectrum  of  ©^(a-X)  showing  the  v’  «.  0,  v"  «  0  band. 

Figure  4.  Low  resolution  spectra  of  0  (a-X)  and  0  (b-X).  (Ref.  22) 

2  2 
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Table  1.  Spectroscopic  constants  for  0^  X,  a,  and  b  states. 
(All  quantities  reported  in  wavenunbers,  cm”') 


State 

T 

e 

u 

e 

Ci>  X 

e  a 

<i)  y 

a  a 

B 

a 

a 

a 

x^- 

0 

1580.36 

12.073 

0.0546 

1 . 44566 

0.01579 

a}t 

7918.1 

1509.0 

12.0 

• 

1.4264 

0.0171 

9 

13195.2 

1432.687 

13.950 

-0.01375 

1 . 40041 

0.01817 

2.2  THE  EPR  SPECTROSCOPY  OF  0  (a)  AND  0,(X) 

- 2 - 2 - 

Electron  Paramagnetic  Resonance  or  Electron  Spin  Resonance  (ESR)  spectro¬ 
scopy  usually  refers  to  microwave  transitions  between  energy  levels  produced 
by  the  application  of  a  magnetic  field  on  an  unpaired  electron  (Refs.  23-24). 
Unfortunately,  the  term  (EPR)  has  been  applied  to  magnetically  induced  transi¬ 
tions  which  are  not  necessarily  paramagnetic  transitions,  and  the  terms  EPR 
and  ESR  have  been  used  interchangeably.  The  simplest  example  of  EPR  spectro¬ 
scopy  is  that  of  a  free  electron  in  the  presence  of  a  uniform  magnetic  field. 
An  electron  possesses  a  magnetic  moment,  fi,  due  to  its  intrinsic  spin,  S. 

This  spin  may  be  oriented  with  or  against  the  applied  field,  H,  yielding  the 
two  spin  states  denoted  by  the  quantum  number,  m^  =  ±  1/2.  The  energy  degen¬ 
eracy  of  these  two  spin  states  is  removed  by  the  magnetic  field.  The  energy 
of  a  magnetic  dipole  in  a  magnetic  field  is  given  by  Equation  1: 

E«-|i«H  =  pH*sgumH  (1) 

Z  B  S 

where: 

■  component  of  the  magnetic  dipole  in  the  direction 
of  the  magnetic  field 

u  ■  Bohr  magneton  ■  9.27  x  lO”*'  ergs/gauss 

B 

g  m  spectroscopic  splitting  factor  »  2.(X)2 
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Therefore,  the  energy  difference  between  the  two  spin  states  is 

AE  «  hw  »  g  u  H  (Ab  ■  1)  (2) 

P  9 

Ihider  the  application  of  a  magnetic  field,  the  electron  will  undergo  a  "spin 
flip"  at  frequency  i>  =  (g  u/h)  H.  The  spectroscopic  splitting  factor  is 

P 

defined  as  g  »  (h/u  )  dw/dH.  Usually  the  derivative  dv/dH  is  a  constaiit.  The 

P 

value  of  gp  /h  for  a  free  electron  is  2.8  MHz/G,  and  thus,  for  magnetic  fields 

P 

on  the  order  of  kllogauss,  the  transition  is  in  the  microwave  (gigahertz) 
region  of  the  electromagnetic  spectrum. 

A  typical  Em  spectrometer  operates  at  a  fixed  mlcromve  frequency,  v,  while 
the  magnetic  field  strength  is  scanned  across  the  resonance  described  by 
Equation  2.  When  the  magnetic  field  satisfies  the  resonance  condition,  the 
microwave  radiation  is  absorbed  by  the  sample  under  study  and  the  transmitted 
intensity  decreases.  Phase  sensitive  detection  is  often  employed  to  improve 
the  detection  limits.  Figure  5  schematically  illustrates  the  microwave  field 
strength  as  a  function  of  time.  If  the  modulation  amplitude  is  small  compared 
to  the  line  width,  then  the  derivative  of  the  absorption  profile  is  observed 
(Ref.  23).  Typical  EPR  signals  for  a  ground  state  oxygen  gas  sample  are  shown 
in  Figure  6. 

This  simple  analysis  of  EPR  spectra  applies  to  gas  phase  molecules  with  a 
single  unpaired  electron.  The  EPR  spectra  can  be  obtained  from  other  systems. 
For  example,  an  oxygen  molecule  has  two  elec'-rons  in  the  unfilled  lx  anti- 
bonding  orbital  and  its  EPR  spectrum  is  more  complex.  Indeed,  for  0^(a  A), 
the  total  electronic  spin  is  zero,  and  thus,  0^(a)  has  no  spin  induced  mag¬ 
netic  moment.  In  addition,  a  single  spin  flip  in  the  ground  state,  0  (^), 

2 

involves  a  chauige  of  electronic  state  (Fig.  2).  A  more  generalized  view  of 
angular  momentum  is  required  to  understand  the  EPR  spectrum  of  molecular 
oxygen. 

The  EPR  spectrum  for  0  (a*A)  has  been  analyzed  In  detail  (Refs.  7,9).  For  the 

a  A  state  of  oxygen,  the  spin  angular  momentum  is  zero  and  for  there  is 

2 
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Magnetic  Field  Amplitude 
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Figure  5.  Phase  sensitive  detection  and  field  modulation  parameters. 
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a)  Observed  ESR  signal  showing  derivative  of  absorption  feature. 


b)  First  integral  of  ESR  signal 
showing  broadened  absorption 
profile. 


c)  Second  integral  of  ESI  signal  whose 
asyntope  provides  area  of  absorption 
feature. 


Figure  6.  Typical  ESR  signals  observed  for  ground  state  molecular  oxygen. 
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no  nuclear  spin.  Thus,  the  total  angular  momentun,  J,  is  simply  the  sun  of 
the  orbital  angular  momentum,  A,  and  the  rotational  angular  momentum,  K,  so 
that  J  "  A  +  The  orbital  angular  momentum  is  related  to  the  magnetic 
moment  by 

M  Mg  A  (3) 

where  the  spectroscopic  splitting  factor  is  g^  =  1.0,  since  only  orbital 
angular  momentum  is  involved.  The  energy  of  the  magnetic  interaction  is 
E  *  -y  •  H,  and  for  AM^  =  1,  the  energy  splitting  is 

AE  »  gj  H  (4) 

2 

where  the  spectroscopic  splitting  factor  g  *  A  /J(J+1)  is  obtained  from  the 

J 

projection  of  the  orbital  angular  momenttim  onto  the  axis  of  the  magnetic  field 
(Ref.  9).  The  lowest  energy  angular  momentum  level  of  0  (^A)  is  A  =  2,  K  =  0 
and  J  *  2  which  yields  a  value  for  g  of  2/3.  For  J  =*  2  there  are  five  values 
for  the  projection  of  J  onto  the  axis  of  the  magnetic  field,  M  =  -2,  -1,  0, 

J 

1,  2.  According  to  Equation  4,  these  states  are  equally  spaced  in  energy. 
However,  in  second  order  perturbation  theory  the  transition  degeneracy  is 
lifted  and  the  energy  splittings  for  AM^  »  1  become  (Ref.  9) 

AE  «  g,  u  H  -  (2M +1)  A^  g,^  u  ^  H*/  378B  (5) 

J  p  J  LB  O 

where  is  the  rotational  spectroscopic  constamt  for  the  vibrational  level 
v'  =  0  providing  the  rotational  energy  E  =  B  J(J+1).  Figure  7  illustrates 

rot  o 

the  splitting  of  these  M  levels  in  the  presence  of  a  magnetic  field  for  both 

J 

the  J  >  2  and  J  »  3  levels  of  O^(a^A).  The  J  >  3  level  has  a  degeneracy  of 
2J4-1  *  7.  The  rotational  energy,  B  «  1.41808  cm~^  (Ref.  22),  is  large  com- 

O 

pared  to  the  magnetic  splitting.  Using  Equation  4,  it  follows  that  for  a 
fixed  field  strength,  the  splitting  between  M  states  in  the  J  «  3  manifold  is 
smaller  than  the  splitting  in  the  J  «  2  manifold.  To  access  the  transitions 
in  the  J  >  3  manifold  using  a  klystron  frequency  on  the  order  of  10  GHz,  a 


12 


13 


WL-TR-90-85 


static  nagnetic  field  of  a  20  KG  is  needed,  which  is  beyond  the  capacity  of 
the  system  used  in  this  study. 

The  magnetic  field  strengths  experimentally  observed  by  A.  M.  Falick  (Refs. 
6,9)  for  the  four  AM  *  1  transitions  of  J  »  2  at  a  microwave  frequency  of 
9.29  GHz  are  reported  in  Table  2.  The  observed  values  for  hv/gyt^H  may  be 
used  to  calculate  the  resonant  field  strengths  at  other  microwave  frequencies 
from  Eqxiatlon  4.  Also  listed  in  Table  2  are  the  predicted  field  strengths  for 
the  microwave  frequency  of  8.917  GHz  used  in  the  present  study. 

A  more  complete  solution  for  the  Hamiltonian  of  a  rotating  molecule  in  an 
external  magnetic  field  may  be  found  in  Reference  7.  The  effect  of  nuclear 
spin  and  a  transformation  of  the  Hamiltonian  to  handle  the  coupling  of  differ¬ 
ent  states  is  presented.  The  line  positions  for  O^(a^A),  J  »  2  smd  J  »  3 
were  observed  with  considerably  greater  resolution  than  in  the  previous  study 
(Ref.  5)  and  yield  a  precise  set  of  rotational  and  magnetic  constants  that  are 
entirely  self  consistent.  The  improved  resolution  provided  by  Reference  7  is 
not  required  for  the  relatively  low  resolution  work  to  be  reported  in  this 
document.  The  relative  Intensities  of  the  J  »  2,  AM  i  transitions  have 

J 

been  calculated  and  are  in  the  ratio  2: 3: 3: 2  for  the  M  »  -2,  -1,  0,  and  -t-l 

J 

initial  states,  respectively  (Ref.  9). 

The  EPR  spectrum  of  ground  state  oxygen  has  also  been  studied  in  detail 
(Refs.  5,25).  For  the  X^Z  ~  electronic  state  of  oxygen,  the  orbital  angular 

9 

momentum  is  A  »  0,  and  the  total  electronic  spin  is  Z  >■  1.  Coupling  the  elec¬ 
tronic  spin  to  the  molecular  rotation,  {C,  produces  the  total  angular  momentum, 
i  *  Z  K.  For  K  2  1,  each  rotational  state  produces  three  values  of  total 
angular  momentum,  J  »  K  1,  K,  and  K  -  1.  The  configurations  leading  to  the 
three  triplet  spin  states  are  shown  in  Figure  2.  Coupling  the  magnetic  moment 
of  0^  to  an  external  magnetic  field  splits  the  J  manifold  into  2J-*-! 
nondegenerate  states. 
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Table  2.  Microwave  trsmsltlons  of  0^(a)  in  a  magnetic  field.  (Ref.  6) 


Transition 

Observed  Field 

at  V  -  9.29 

(Gauss) 

GHz 

hv/p  H 

B 

(observed) 

Predicted  Field  (Gauss) 

at  V  e  8.917  GHz 

10,090.3 

0.65596 

9,712.6 

9,988.7 

0.66264 

9,614.7 

9,885.6 

0.66956 

9,515.4 

B9QI 

9,781.4 

0.67640 

9,419.1 

The  energy  of  the  magnetically  split  K  »  3  rotational  level  as  a  function  of 
magnetic  field  strength  is  shown  in  Figure  8.  The  J  -  K  ±  1  lines  are  split 
from  the  J  »  K  lines  due  to  spin-spin  interactions  of  about  60  GHz  (Ref.  5). 
Spin-orbit  and  spin-rotation  Interactions  are  smaller  and  split  the  J  ^  K  +  1 
and  J  s  K  -  1  lines  by  about  -4  to  7  GHz,  depending  on  the  rotational  state,  K 
(Ref.  S).  A  list  of  observed  field  strengths  for  various  transitions  at  a 
microwave  frequency  of  9.000  GHz  are  provided  in  Table  3.  Also  listed  are  the 
relative  transition  strengths  and  values  for  dv/dH  »  gp^/h  used  to  calculate 
the  resonant  field  strength,  H,  for  any  microwave  frequency,  v.  Due  to  spin- 
rotational  coupling  and  the  large  interaction  between  unpaired  electrons  in 
ground  state  oxygen,  the  resonant  frequency  is  not  linearly  dependent  on  0»g- 
netic  field.  That  is,  dv/dH  is  not  a  constant.  In  addition,  the  transition 
probabilities  vary  with  the  freq[uency.  These  parameters  have  been  measured  as 
a  function  of  frequency  in  the  range  8.5  to  9.6  GHz  (Refs.  5,25).  The  line 
positions  and  strengths  can  be  calculated  to  within  ±2.5  percent  (Ref.  5). 

The  change  in  dw/dH  and  the  transitions  probabilities,  with  frequency 

are  also  given  in  Table  3. 

2.3  RELATIONSHIP  OF  OBSERVED  EPR  SIGNALS  TO  CONCENTRATIONS 


Ground  state  oxygen  has  often  been  used  as  a  standard  to  obtain  absolute 

concentrations  from  EPR  spectra  (Refs.  6,26).  Absolute  02(a*A)  concentration 

are  obtained  by  observing  the  relative  EPR  signals  for  0  (^A)  and  0  (^)  in  a 

2  2 

pure  oxygen  flow.  Given  the  total  gas  pressure,  temperature  and  the  ratio. 
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Table  3.  Frequency  dependence  of  transition  probabilities  for  the  principal 

0,  lines.*  (Ref.  5) 

2 


transition 

A (dv/dB) 

A4[S  IVGHz 

X 

K 

J 

M 

8.“  0« 

dv/dB® 

4[S  f 

X 

e  a/ 

0 

e  Av** 

0 

5 

4  -» 

6 

1  -»  2 

5094 

1.817 

0.3674 

0.045 

0.0666 

1 

1 

1  -»  0 

5343 

1.945 

0.7354 

0.044 

0.0241 

1 

2 

1^2 

5814 

1.727 

1.2010 

0.033 

0.0242 

1 

2 

0  ->  1 

6393 

1.429 

1.4862 

0.000 

-0.0021 

5 

6  -> 

4 

-2  -»  -1 

6237 

1.730 

0.6265 

0.045 

0.0725 

3 

4 

-1^0 

8293 

1.660 

1.4868 

0.076 

0.0238 

1 

2 

-1^0 

6910 

1.297 

1.2633 

0.011 

-0.0091 

3 

4 

3^4 

9242 

1.294 

0.9217 

0.060 

0.0431 

*Calculation  valid  to  1  1  in  last  figure  for  resonant  frequency  from  8.5  to 
9.6  GHz 

**l  0«  »  O.l  mT 

*MHz/0#  (1  MHz/O#  -0.1  GHz/T) 

Changes  in  dv/dB  and  4(S^]^  per  GHz  in  resonant  frequency 


r  ■  lO  (*A)]/[0  (^)]  and  assuming  the  presence  of  only  0  (a)  and  0  (X)  yields 
•  *  0  2  2 

the  relationships: 


IN^I  »  P/kJ  -  lO^CX)]  ♦  (O^Ca)! 


(1+r)  (02(X)] 

(7) 

r/(l+r)  (P/k  T) 

B 

(8) 

The  extremely  small  errors  associated  with  assuming  only  0  (a)  and  0  (X)  are 

2  2 

present  in  the  flow  will  be  addressed  in  Section  3.3. 


Berg,  J.O. ,  "EPR  Calibration  Report",  Chemical  Oxygen-Iodine  Laser 
Advanced  Technology  (COIL)  USAF  Contract  No  F29601-80-C-0028,  12  April  1982. 
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The  absorbed  microwave  power  is  related  to  the  number  density  of  the  observed 
state,  with  the  proportionality  constant  depending  on  the  trauisltion  probabil¬ 
ity  and  the  statistical  population  distributions;  specifically,  the  concentra¬ 
tion  of  state  i,  [NJ,  is  given  by  (Ref.  23): 

(Njl  =  (2k^T/hv)  expC-E^/k^T)  (9) 

where: 


=  energy  of  state  i 

Z  B  molecular  partition  function 
Ip  I  a  transition  probability  from  state  i  to  state  J 

-  Instrument  constant.  Independent  of  species  or  state 
I  »  (1/G  H  P  j"a'dH'dH  =  integrated  EPR  signal  for  state  i 

1  lapOOl 

(H)  ■  dot/dH  »  observed  EPR  signal  for  state  i 
a^(H)  •  microwave  absorption  coefficient  for  state  1 
Gj  »  EPR  gain  setting  for  transition  (i,J) 
s  EPR  magnetic  field  modulation  amplitude 
P^  *  incident  microwave  power 


Thus,  the  ratio  r  *  (0^(a)]/[02(X)] ,  is  simply  the  ratio  of  the  observed 
intensities  multiplied  by  a  relative  intensity  factor,  (x)'^^o  {  )’  ^ 

relative  gain  setting,  G  .  /G  ,  22 


I0,(a)]/0,(X)] 
2  2 


)  (S  /S  )  (G  /G  ) 
0  yc)  0  (•)  0  (X^  0  (m)^ 


(10) 


The  relative  intensity  factor,  S  ,„,/S  ,  is  a  collection  of  state  depen- 

2  2 

dent  constants  from  Equation  9.  Based  on  the  reported  transition  probabili¬ 


ties,  the  relative  Intensity  factor  for  the  0^(a)  J  ■  2,  ■  0  -»  1  transition 
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with  respect  to  the  0  (X)  K«3,  J-4.  M  b3-»4  transition  at  a  microwave 

2  J 

frequency  of  9.347  GHz  and  temperature  T  »  300  K  has  been  determined  as 

(S  /S  ,  °  0.624  (Ref.  4).  This  value  is  approximately  linearly 

o_cx)  o 
2  2 

dependent  on  temperature  due  to  the  partition  function  (Ref.  5).  Due  to  the 

nonlinear  dependence  of  the  0^^^^  transition  probabilities  on  frequency,  this 

value  also  depends  slightly  on  frequency.  This  small  frequency  effect  can  be 

accounted  for  with  the  values  of  A|p^^|/Aw  provided  in  Table  3.  The  value 

(S  /S  )  =  0.624  is  use..  '  sgar  Jless  of  microwave  frequency. 

2  2 

The  broadening  of  the  absorption  profile,  a(H),  can  have  significant  effects 
on  the  double  integration  require i  to  obtain  the  intensities,  I^.  Several 
mechanisms  including  collisional  relaxation  and  modulation  will  broaden  the 
absorption  feature.  Often,  observed  profiles  exhibit  both  homogeneous  and 
inhomogeni 3US  behavior.  Volght  profiles  may  be  used  for  these  convolved  cases 
(Ref.  27).  Lorentzian  line  shapes  (characteristic  of  homogeneous  broadening) 
decay  slower  than  Gaussian  line  shapes  from  resonance,  (indicative  of  inhomo¬ 
geneous  broadening).  A  significant  error  can  be  obtained  when  tnincating  the 
integration  of  a  Lorentzian  profile  (Ref.  23).  However,  correction  factors 
for  this  trxincation  can  be  easily  obtained  (Ref.  23).  The  effects  of  modula¬ 
tion  broadening  and  saturation  have  been  analyzed  (Refs.  28-32). 

2.4  TRANSPORT  PROPERTIES  OF  0  (a) 

— '  —  ■  2 - 

The  design  of  the  calibration  setup  did  not  permit  the  microwave  absorption  to 
be  measured  at  the  same  location  as  the  infrared  emission  detection.  Thus  the 
effects  of  deactivation  during  transport  of  the  singlet  oxygen  must  be  con¬ 
sidered.  Equation  22  provides  the  basis  for  this  correction  to  the  detector 
calibration,  as  described  in  Section  4.0.  The  efficient  transport  of  singlet 
oxygen  is  made  possible  by  both  its  long  radiative  lifetime  and  extremely  slow 
collisional  deactivation  rates.  The  gas  phase  kinetics  of  singlet  oxygen  are 
well  known  (Ref.  33).  There  are  several  paths  for  energy  removal  including 
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radiation,  homogeneous  gas  phase  quenching,  wall  deactivation,  and  energy 
pooling: 


0,(*A) 

2 

^  0,(\) 
2 

+ 

hv 

(11) 

0  (*A)  +  M 

2 

-»  0,(^) 
2 

+ 

M 

(12) 

0,('a)  +  0,('a) 

2  2 

-»  0,(^Z) 

+ 

o,(^) 

(13) 

0,(‘a)  +  0,('a) 

2  2 

+ 

o,(^) 

2 

(14) 

0,('a)  +  0,(^A) 

2  2 

+ 

0  (^Z) 

2 

(15) 

0  (‘a)  +  wall 

2 

(16) 

In  addition,  the  deactivation  of  0^(^£)  can  remove  energy  from  the  system 
through  the  following  reactions: 


0^(‘Z) 

+  hv" 

(17) 

+  M 

+  M 

(18) 

+  M 

(19) 

0  (^Z)  +  wall 

2 

^  Ojj(’i) 

+  M 

(20) 

+  M 

(21) 

Deactivation  by  radiative  paths  may  be  neglected  due  to  the  extremely  long 

lifetimes  of  0^(a)  and  0^(b).  Recently,  Knickelbeln  and  coworkers  (Ref.  34) 

confirmed  the  long-standing  supposition  that  0  (^Z)  is  quenched  to  0  (^A) 

3  2  2 
rather  than  to  0^(  Z).  Thus,  reactions  19  and  21  may  be  neglected.  An 

analytical  solution  for  the  concentration  of  singlet  oxygen  as  a  function  of 

time  can  be  obtained  when  the  O^C^Z)  concentration  has  achieved  steady-state 

conditions  (Ref.  35): 

[0,(‘a)]"‘  «  {  I0,(^A))  }  exp(Bt)  -  A/^  (22) 

2  2  O 
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where: 

[0  (^A)]  *  concentration  of  0  (^A) 

2  2 

[0  (*A)]  =  concentration  of  0  (  A)  at  t  «  0 

2  0  2 

A  .  »  N.  *  ''is 

10^(*Z)1  =  k,3  tO^C'A))*  /  (k33  ♦  k^^lW) 

[M]  =  concentration  of  arbitrary  quencher;  predominantly 
water  and  0^(X) 

The  steady-state  condition  for  singlet-sigma  oxygen  is  attained  rapidly.  Note 
that  the  relative  importance  of  quenching  (Equation  12}  and  pooling  (Equation 
13)  depends  critically  on  the  excited  oxygen  yield.  In  typical  flow  tube 
experiments  quenchiiig  dominates  the  deactivation,  whereas  pooling  dominates 
under  laser  device  conditions. 

In  the  development  of  chemical  generators  of  singlet  oxygen  a  correlation 
between  singlet  oxygen  concentration  and  total  pressure,  partial  pressure  of 
singlet  oxygen,  temperature  and  total  flow  rate  has  been  developed  amd  termed 
the  "chi  correlation"  (Ref.  12).  This  relationship  Includes  only  second  order 
kinetic  deactivation  and  should  not  be  applied  to  microwave  sources  of  singlet 
oxygen  where  the  yield  is  low  and  ground  state  oxygen  quenching  of  0^(a)  is 
important. 
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3.0  EXPERIMENTAL  APPROACH 

3.1  CALIBRATION  APPARATUS 

A  schematic  of  the  experimental  apparatus  for  calibrating  the  infrared  emis¬ 
sion  detector  using  EPR  spectroscopy  is  shown  in  Figure  9.  The  apparatus  con¬ 
sisted  of  four  major  subsystems:  (1)  a  flow  tube  and  source  of  singlet  oxygen, 
(2)  an  optical  collection  system  and  infrared  detector,  (3)  an  EPR  spectro¬ 
meter,  and  (4)  a  data  acquisition  system.  Each  of  these  subsystems  will  be 
described  in  detail,  followed  by  a  discussion  of  important  design  Issues. 

A  moderate  gas  velocity  (900-4500  cai/s)  flow  tube  was  designed  to  generate  and 
transport  a  laminar  flow  of  singlet  and  ground  state  oxygen  to  a  calibration 
volume  where  the  emission  from  0^(a)  was  detected  and  to  the  cavity  of  an  EPR 
spectrometer  where  the  relative  0^(a)  and  0^(X)  concentrations  were  observed. 

A  Kinney  model  850  mechanical  vacuum  pump  (800  CFM)  assisted  by  a  Roots  model 
MF2700  blower  (24(X)  G^)  was  used  to  establish  a  flow  rate  of  50  to  1000  seem 
of  pure  oxygen  controlled  by  a  Unit  Instruments  model  URS  100-5  mass  flow  con¬ 
troller  with  a  model  UFC  1000  sensor  operating  in  the  range  of  0  to  2000  seem. 
Pressures  as  low  as  17  mtorr  and  leak  rates  <  40  mtorr/min  were  achieved. 
Pressures  in  the  range  of  50  to  1500  mtorr  were  recorded  with  MKS  model  390HA 
1-torr  and  model  122AA  10-torr  capacitance  manometers  calibrated  for  zero, 
linearity  and  span  on  a  MKS  model  PVS-2A-0100  with  270B  signal  conditioner 
calibration  stand.  A  linear  regression  analysis  provided  a  slope  of  1.0020 
and  an  intercept  of  0.002.  Pressure  readings  were  accurate  to  within  0.5 
percent.  A  manifold  was  constructed  to  sample  the  pressure  at  each  of  eight 
different  locations  with  the  same  sensor.  The  zero  calibration  of  the  pres¬ 
sure  transducer  was  checked  regularly  on  a  2-in  oil  diffusion  pump. 

Total  flow  rates  were  measured  and  actively  controlled  by  the  controller 
operating  in  the  range  of  0  to  2000  seem.  The  manufacturer’s  calibration  for 
nitrogen  gas  was  corrected  for  oxygen  flow  by  the  square  root  of  the  ratio  of 
heat  capacities,  ^ "  0.996.  The  gas  temperature  was  measured 
with  Medtherm  type  T  gas  phase  thermocouples  to  an  accuracy  of  ±  0.  l"c.  The 
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temperature  did  not  depend  on  flow  conditions  or  microwave  power  and  the  room 
temperature  was  maintained  at  a  constant  21  **0.  The  measured  flowing  gas  tem¬ 
perature  was  always  T  »  21  ±  1  **0.  The  locations  for  pressure,  temperature 
and  mass  flow  measurements  are  Indicated  In  Figure  10. 

Ultra-high  purity  oxygen  was  passed  through  a  Opthos  Instruments  gas  nitrogen 
cooled  Evanson  microwave  cavity  via  1.27-cm  diam  quartz  tubing  to  generate  5 
to  14  percent  O^(a^A).  The  cavity  was  excited  by  a  Klva  model  MPG-4  2480-  MHz 
magnetron  generator  with  70  W  forward  and  1  to  2  W  reflected  power.  A  Wood’s 
horn  situated  Just  downstream  of  the  microwave  cavity  was  used  to  reduce  the 
scattering  of  the  microwave  generated  emission. 

The  reaction  0  +  NO  N0_  hv  leads  to  broad  band  (visible  and  IR)  emission. 

2 

The  emission  Is  within  the  1.27  pm  band  pass  of  the  detector  and  can  cause 
calibration  errors.  The  microwave  cavity  produces  both  oxygen  and  nitrogen 
atoms.  To  eliminate  oxygen  atoms,  a  mercury  oxide  coating  was  applied  to  the 
walls  of  the  microwave  exit  channel  to  Induce  oxygen  atom  recombination.  This 
method  has  been  used  by  several  research  groups  (Refs.  37,  39).  Air  leaks 
upstream  of  the  microwave  cavity  are  the  only  source  of  nitrogen.  There  were 
no  air  leaks  detected  upstream  of  the  microwave  discharge  with  the  leak  detec¬ 
tor  used  In  this  work.  Consequently,  no  recombination  emission  was  observed 
under  visual  inspection.  A  more  detailed  analysis  of  this  emission  intensity 
Is  provided  In  Section  3.3. 

The  oxygen  was  transported  to  the  calibration  volume  through  2.21-cm  Inner 
diameter  Pyrex  tubing.  One  of  three  optical  detection  blocks  (Fig.  11)  was 
connected  to  the  Pyrex  tubing  with  Nacom  2.S4-cm  diam  Teflon  fittings  to 
monitor  the  concentration  of  02(a^A)  delivered  by  the  microwave  generator. 

TWO  additional  detection  blocks,  B1  and  BZ,  upstream  and  downstream  of  the  EPR 
cavity  were  used  to  monitor  singlet  oxygen  deactivation.  To  maintain  laminar 
flow  and  prevent  axial  pressure  gradients,  entrance  ramps  with  an  expansion 


*  Ibid,  pg  17. 
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Figure  11.  Drawing  of  optical  detection  block. 
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half-angle  of  12^  were  designed  to  transition  the  flow  from  the  1-in  dlaa 
tubing  to  the  17.8-  by  4.0-cn  cross  section  of  the  calibration  volume.  The 
interior  of  the  entrance  ramps  and  calibration  duct  were  coated  with  identical 
black  Kynar  material  used  in  the  diagnostic  duct  of  the  COIL  device  (this 
material  is  used  on  COIL  to  reduce  wall  deactivation  of  singlet  oxygen). 

Thus,  any  differences  in  light  scattering  between  the  calibration  duct  and 
the  COIL  device  were  eliminated. 

It  was  possible  to  detect  the  A  *  1.268  pm  emission  from  0^(aA,  v'  *  0)  to 

0  (X^,  s  0)  at  any  of  the  three  detector  blocks  (Bl,  B2,  B3)  or  any  of 
2 

the  three  windows  in  the  viewing  volume  hardware.  The  optical  collection 
system  included  an  Oriel  quartz  f  «  75  mm  collimating  lens  with  a  limiting 
aperture  of  1.07  cm.  The  emission  from  the  collimated  viewing  volume  was 
focused  onto  the  end  of  an  Oriel  3-mm  diam  by  9.1-m  long  glass  fiber  optic 
bundle  for  transmission  to  the  detector  assembly.  According  to  the  manufac¬ 
turer  the  fiber  optic  introduced  an  optical  loss  of  ~  3  percent  per  meter  at 
A  «  1.268  pm.  To  reduce  the  sensitivity  of  the  detection  system  to  changes  in 
collection  efficiency,  transmission,  optical  focus,  and  geometry,  the  complete 
optical  and  detector  system  was  designed  as  a  single  unit,  as  shown  in  Figure 
12.  An  interference  filter  centered  at  A  «  1.268  pm  with  a  full-width  at 
half  maximum  of  10  nm  was  used  to  limit  the  the  detected  emission  to  the 

single  vibrational  transition  of  0  (a*A,  v'  «  0)  to  0  (X^,  v"  »  0).  The 

2  2 

average  trauismission  of  the  filter  was  26.4  percent  for  the  y'  »  0,  v*  «  0 
transition.  The  detector  housing  could  accommodate  neutral  density  filters  to 
attenuate  signals  incident  at  the  detector.  The  optic  probe  movint  included  a 
shutter  to  eliminate  all  light  incident  on  the  detector  when  necessary.  The 
output  from  the  fiber  optic  bundle  was  focused  with  a  f  »  75  nm  lens  to  over¬ 
fill  the  active  area  (25  mm^)  of  the  IR  detector. 

A  liquid  nitrogen  (LN  }  cooled  intrinsic  germanium  detector  manufactured  by 
the  Applied  Detector  Corporation,  model  403L,  was  used  to  detect  the  Infrared 
emission.  The  radiant  flvix  incident  on  the  detector  due  to  the  0  (a)  emission 

14  * 

for  a  concentration  of  10  molecules  per  cubic  centimeter  was  calculated  to 
be  m  10~*^  U.  This  detector  has  a  large  responsivlty,  R  *  7  x  10*  V/W,  and 
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Figure  12.  Scheiaatlc  diagram  of  0^(a)  optical  detection  system. 
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moderate  noise  equivalent  power,  NEP  a  8.6  x  10  W,  so  that  the  detectivity 

4  1/2  2 
is  large,  D  «  5  x  10  m-Hz  /W.  The  active  area  of  the  detector  is  25  am 

and  has  a  linear  response  in  the  chopper  frequency  range  10  to  500  Hz.  The 

responsivity  peaks  at  1.6  pm  and  drops  by  «  60  percent  at  1.268  jim. 

The  detector  required  a  power  supply  (ADC  model  PS-1  or  PS-2)  operating  at 
-250  VDC.  The  detector  DC  output  voltage  served  two  functions:  indicating 
proper  operating  temperature  and  optical  saturation  of  the  detector.  When 
cooled  to  77  K  the  detector  output  DC  voltage  read  -1.25  V.  A  DC  voltage  near 
9  V  indicated  optical  saturation.  The  detector  signal  due  to  background  light 
was  vindetectable  with  a  dark  signal  of  a  50  pV.  Further  characterization  of 
the  detector  is  provided  in  Section  3.3. 

A  single  IR  detector  was  used  to  record  up  to  three  different  emission  signals 
by  chopping  the  light  at  three  different  frequencies  and  using  a  trifurcated 
fiber  optical  bundle  to  couple  the  emission  to  the  detector's  active  area. 
Phase  sensitive  detection  was  used  to  increase  the  detection  limit  and  signal 
to  noise  ratio.  Two  PAR  model  192  dual -frequency  mechanical  choppers  provided 
the  distinct  frequencies  of  94.33,  79.39,  and  55.5b  Hz.  The  detector's  AC 
output  voltage  was  processed  by  Ithaco  model  395  Heterodyne  Narrowband  Volt¬ 
meters  with  auto-phasing  capability.  Typical  signals  ranged  from  0.3  to  2  mV 
for  a  time  constant  of  0.3  s. 

A  Varlaui  E-Llne  model  E-112  Electron  Paramagnetic  Resonance  Spectrometer 

(Ref.  40)  was  used  to  measure  the  relative  0  (a)  and  0  (X)  concentrations.  A 

2  2 

Varian  V-7300-12  in  electromagnet  established  a  variable  magnetic  field 
strength  of  0  to  10,000  G.  The  magnet  was  driven  by  a  V-7700  5-to  90-A  13-kW 
power  supply.  The  magnetic  field  was  modulated  at  100  kHz  with  a  peak-to-peak 
amplitude  of  0.5  to  4.0  G  by  an  E-207  High  Frequency  module.  A  model  E-102 
Microwave  Bridge  with  Klystron  was  the  source  of  the  X-band  microwave  radia¬ 
tion  with  a  frequency  of  8.8  to  9.6  GHz  and  a  maximum  power  of  200  mU.  A 
rectangular  cross  section  waveguide  with  a  60-dB  variable  attenuator  was  used 
to  transport  the  microwave  radiation  to  the  model  E-231  cavity.  The  rectangu¬ 
lar  TE  cavity  provided  an  unloaded  quality  of  Q  ■  7000.  Typically,  the  EPR 
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spectrometer  was  operated  at  the  fixed  microwave  frequency  8.92  GHz  while  the 
magnetic  field  strength  was  scanned  by  20  or  40  G  over  either  the  O^CX^) 
K>3.  J«4,  M  »3-»4  transition  at  9171  G  or  the  0  (a^A)  J«2,  M  »0-»l 

J  2  J 

transition  at  9515  G.  A  gain  of  250  for  0  (X)  and  of  3200  for  0  (a)  were 

2  2 

typically  used  to  provide  measured  signal  to  noise  ratios  of  500:1  and  100:1, 

respectively.  The  detected  EPR  signal,  the  derivative  of  the  absorption  as  a 

fimctlon  of  magnetic  field  a’ (H).  was  recorded  on  an  E-200A  oscilloscope, 

E-080A  X-Y  plotter,  and  the  computer  data  acquisition  system  A  typical  40  G 

scan  required  2  min  for  0  (^A)  and  1  min  for  0  (^). 

2  2 

The  calibration  data  were  acquired  and  stored  on  a  Zenith  Z-248  personal 

computer  system  for  further  data  analysis.  Pressure,  temperature,  mass  flow 

rate,  IR  detector  signals  at  each  of  three  locations  (B2,  Calibration  Volume, 

and  B3).  and  a* (H)  were  recorded  using  the  ASYST  laboratory  software  package. 

Data  were  acquired  at  a  rate  of  1000  points  per  scan.  A  series  of  data  files 

for  various  pressures,  flow  rates,  0  (a)  microwave  powers,  and  EPR  spectro- 

2 

meter  settings  were  obtained  for  each  detector  calibration.  The  first  and 
second  Integrals  of  the  observed  EPR  signal,  a'(H),  were  computed  and  dis¬ 
played.  The  areas  were  obtained  using  a  Simpson  Integrating  routine  after 
subtraction  of  the  baseline.  The  baseline  was  calculated  as  an  average  of 
all  the  a'  (H)  data  points,  as  suggested  in  Reference  23. 

3.2  CALIBRATION  PHILOSOPHY 

To  improve  the  calibration’s  validity,  every  attempt  was  made  to  make  the 
calibration  conditions  as  similar  as  possible  to  those  expected  on  the  COIL 
where  the  detector  was  to  be  used.  The  detector  was  calibrated  on  a  flow 
apparatus  with  identical  viewing  volume,  reflective  surfaces,  viewing  windows, 
etc.  as  the  COIL  diagnostic  duct.  This  eliminated  potential  problems  of 
viewing  volume  differences  between  the  calibration  duct  and  the  COIL  duct. 

For  large  viewing  volumes,  it  may  be  impractical  to  calibrate  the  detector  on 
an  exact  replica  of  the  laser  hardware.  Large  volumes  imply  long  transport 
times  and  greater  deactivation  of  the  singlet  oxygen.  As  a  result,  the  signal 
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to  noise  ratio  of  the  O^Ca)  EPR  signal  was  substantially  reduced  and  the 
required  correction  for  deactivation  is  large  and  uncertain.  Indeed,  these 
errors  exceed  those  caused  by  transferring  the  calibration  to  different  view¬ 
ing  volumes.  The  error  analysis  presented  in  Section  5.0  provides  the  basis 
for  assessing  the  relative  iaportance  of  these  various  sources  of  error  and 
selecting  a  best  calibration  approach  for  a  given  situation. 

A  detector  assembly  was  designed  and  fabricated  as  a  single  unit  to  reduce  the 
sensitivity  of  the  detection  system  to  changes  in  collection  efficiency, 
transmission,  optical  focus,  and  apparatus  geometry.  Once  assembled,  the 
detection  system  remained  intact. 

A  series  of  at  least  10  calibrations  was  performed  on  five  IR  detectors  using 
the  apparatus  and  procedures  described  in  this  document  during  the  period  from 
September  1985  to  December  1988.  This  extensive  set  of  data  provided  a  unique 
opportunity  to  assess  the  reproducibility  and  long  term  stability  of  the  cali¬ 
brations.  A  complete  history  of  the  calibrations  and  data  for  the  application 
of  this  diagnostic  on  a  high  power  CX3IL  device  will  be  forthcoming.*  The 
importance  of  multiple  calibrations  and  regular  detector  sensitivity  tests  to 
assess  changes  in  diagnostic  performance  with  tine  cannot  be  overemphasized. 

3.3  CHARACTERIZATION  OF  CALIBRATION  SYSTEM 

A  detailed  assessment  of  the  operating  characteristics  of  the  calibration 
apparatus  is  absolutely  necessary  to  assess  systematic  and  statistical  errors 
associated  with  the  final  detector  calibration  factor.  A  characterization  of 
the  four  subsystems:  (1)  flow  tube  and  microwave  source  of  singlet  oxygen, 

(2)  optical  collection  system  and  infrared  detector,  (3)  EPR  spectrometer,  and 
(4)  a  data  acquisition  system  is  provided  in  the  following  paragraphs.  With 
careful  experimental  procedures,  the  system  can  be  characterized  with  great 
accuracy. 


To  be  published  as  a  Phillips  Laboratory  technical  report. 
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3.3.1  Flow  Tube  and  0  (a)  Source 
- 2 - 

The  fully  developed,  laminar  flow  in  a  circular  tube  of  constant  area  is 
described  as  Hagen-Pousille  flow  (Ref.  41).  The  mean  velocity  for  such  a 
flow  is 

V  »  -(l/8p)  r^  dP/dx  (23) 

O 

where: 

r  =  the  radius  of  the  tube 

o 

p  -  the  gas  viscosity 

dP/dx  =  the  pressure  gradient  in  the  flow  direction 
From  mass  continuity,  m  -  pvA,  the  total  flow  rate  is 

m  »  -(PnrV8pk  T)  dP/dx  (24) 

o  B 

If  the  driving  force  AP  is  proportional  to  the  head  pressure,  then  the  mass 
flow  rate  will  depend  quadratically  on  pressure.  Figure  13  Illustrates  the 
observed  flow  rates  for  pressures  of  0.1  to  10  tcrr  for  blocks  1  amd  2.  Fig¬ 
ure  13b  shows  data  for  pressure  up  to  1  torr,  the  range  used  during  the  cali¬ 
brations.  Indeed,  the  plots  of  flow  versus  pressure  squared  are  extremely 
linear  in  this  pressure  range  with  correlations  of  r  =  0.9992  and  0.9997  for 
blocks  1  and  2,  respectively.  The  slopes  for  a  given  block  are  different  for 
the  low  and  high  pressure  ranges.  For  example,  the  slope  in  the  low  pressure 
region  for  block  1  is  673.0  sccm/(torr)  as  compared  to  a  slope  of  594.1  seem/ 
(torr)  for  the  full  pressure  range,  about  12  percent  lower.  Similarly,  the 
slopes  in  the  low  and  high  pressure  regions  for  block  2  are  1126.4  seem/ 
(torr)*  and  937.1  sccm/(torr)*,  respectively,  with  the  latter  being  17  percent 
smaller.  The  differences  in  slope  for  a  given  block  are  possibly  due  to  a 
nonlinear  response  of  the  vacuum  station  or  hysteresis  in  the  mass  flow  con¬ 
troller.  A  small,  positive  intercept  in  flow  rate  equal  to  44.6  ±0.1  seem 
is  observed  at  both  blocks.  This  intercept  is  attributed  to  a  zero  offset  in 
the  mass  flow  controller,  with  further  Justification  provided  below.  Due  to 
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a)  High  pressure 


•  BkMk  1  a  Biook  t 


b)  Low  pressure 

Figure  13.  Pusp  curve:  flow  rate  as  a  function  of  pressure 
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frictional  losses,  the  pressure  decreases  as  a  function  of  distance  down  the 
flow  tube  according  to 

P*  -  =  (16pk  T/irr*  )  (m/m)  (x  -x  )  (25) 

2  1  Bo  12 

This  equation  is  obtained  by  solving  Equation  24  for  PdP  and  integrating  over 
a  distance  x^  -  x^  for  a  constant  temperature  T.  Figure  14  demonstrates  that 
the  observed  pressure  drop  between  the  two  optical  detection  blocks  B1  and  B2 
obeys  Equation  25.  Figures  13  and  14  provide  a  good  system  check  for  the  cal¬ 
ibration  of  the  capacitance  manometers  and  flow  meters.  The  fitted  intercept 
of  Figure  13  agrees  favorably  with  that  of  Figure  14.  The  predicted  slope  of 
Figure  14  from  Equation  25,  1656  ±  59  sccm/torr^,  agrees  favorably  with  the 
measured  value,  1620  ±  7  sccm/torr^.  The  uncertainty  in  the  predicted  slope 
arises  primarily  from  the  error  in  determining  the  tube  radius,  r^.  An 
absolute  error  of  <0.1  mm  provides  a  3  percent  uncertainty  due  to  the  fourth 
power  dependence.  Thus,  any  systematic  error  in  pressure  or  flow  rate  other 
than  the  mass  flow  offset  already  identified  must  be  <  4  percent,  and  is  more 
likely  <  1  percent.  It  is  also  demonstrated  that  Equation  25  can  be  used  to 
extrapolate  between  the  two  pressure  readings  at  B1  and  B2  to  obtain  the 
pressure  in  the  EPR  cavity.  The  EPR  cavity  pressure  is  required  to  evaluate 

the  relative  0  (a)  to  0  (X)  concentration  from  Equation  8. 

2  2 

Since  the  flow  rate  is  quadratically  dependent  on  pressure,  the  velocity  must 
be  linearly  dependent  on  pressure.  The  Mach  number  M  *  v/a,  where  a  *  465.5 
m/s  is  the  speed  of  sound  in  oxygen  at  T  »  300  K,  is  linearly  dependent  on 
pressure  as  shown  in  the  plot  of  Figure  15.  The  data  have  been  corrected  for 
a  mass  flow  controller  offset  of  44.6  seem  and  that  successfully  brings  the 
data  to  a  zero  intercept.  The  Mach  number  is  low,  M  s  0.10,  and  the  fluid  is 
incompressible.  This  characterization  of  the  flow  velocity  will  be  required 
to  calculate  transport  times  and  correct  for  the  deactivation  of  singlet 
oxygen,  as  described  in  Section  4.0. 

Characterization  of  the  flow  in  the  calibration  volume  is  complicated  by  the 
variation  of  cross-sectional  area  with  flow  distance.  The  area  profile  for 
the  flow  tube  is  shown  in  Figure  10.  The  observed  pressures  in  the  diagnostic 
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duct  as  well  as  the  pressures  at  Bl,  BZ,  and  B3  for  various  flow  rates  are 
shown  in  Figure  16.  The  predicted  pressures  based  on  a  one  diaensional  aero- 
kinetic  model,  DEACT  (Ref.  42),  are  also  shown  in  Figure  16.  The  observed 
pressures  agree  favorably  with  the  computational  predictions  and  this  computer 
model  will  be  used  in  Section  4.0  to  assess  the  extent  of  singlet  oxygen 
deactivation. 

The  velocity  changes  with  the  cross-sectional  area  to  conserve  mass  flow 
according  to  Equation  26, 

V(x)  »  m/pA(x)  =  m  (k  T/P)  /  A(x)  (26) 

B 

Thus,  the  variation  of  velocity  with  flow  distance,  x,  can  be  obtained  for  any 
flow  condition  from  the  data  presented  in  Figures  10  and  16. 

A  check  for  any  radial  variation  in  pressure  across  the  flow  in  the  diagnostic 
duct  was  performed  by  measuring  the  pressure  at  the  center  of  the  flow  and  at 
the  side  wall.  No  difference  in  pressure  was  observed  over  a  wide  range  in 
flow  rates.  The  sensitivity  of  this  test  was  better  than  0.2  percent.  Thus, 
no  problems  with  nonuniform  flow  fields  or  "Jet-streaming”  were  encountered. 

The  efficiency  of  the  microwave  generator  at  producing  O^('a)  decreased  at 
increased  cavity  pressures,  as  shown  in  Figure  17.  These  data  were  acquired 
over  several  different  test  days.  Absolute  concentrations  were  obtained  from 
the  detector  calibration  reported  in  Section  4.0.  The  scatter  in  the  data  is 
relatively  large,  ±  15  percent.  This  scatter  is  due  to  both  the  lack  of 
reproducibility  of  0^(a)  EPR  absorption  areas  and  the  sensitivity  of  the  EPR 
microwave  generator  to  cavity  tuning.  Also,  no  attempt  was  made  during  the 
calibration  to  reproduce  a  given  yield  (as  might  be  done  by  reproducing  flow 
and  pressure  conditions  and  setting  the  forward  and  reverse  powers  of  the 
microwave  power  supply). 

Also  shown  in  Figure  17  are  the  02(a)  generation  efficiencies  observed  by 
other  researchers.  It  is  quite  typical  to  obtain  5  to  15  percent  of  02(a)  in 
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Figure  16.  Pressure  profile  in  the  flow  direction  within  the  calibration 
duct.  (•)  observed  data,  ( - )  DEACT  model  predictions. 


Figure  17.  Efficiency  of  the  microwave  generator  at  producing  singlet 
oxygen.  (•)  data  from  the  current  study,  (o)  data  from 
Reference  8. 
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such  microwave  systems  (Refs.  4,8).  Such  a  limitation  provides  a  quick  check 

on  the  magnitude  of  any  calibration.  Note  that  the  partial  pressure  of  0  (a) 

2 

continuously  increases  with  pressure  at  least  to  1.5  torr  total  0  pressure. 

2 

The  microwave  source  for  O^(a^A)  is  extremely  stable.  The  IR  signal  was 
observed  by  the  Wang  detector  over  a  period  of  16  min.  No  long  term  change 
in  the  signal  is  observed  and  noise  fluctuations  are  <0.5  percent.  This 
result  implies  that  the  both  the  0^(*A)  and  O^(^)  flows  remain  constant 
during  the  time  required  to  obtain  their  EPR  spectra.  Althoiigh  the  stability 
of  a  standard  microwave  generator  may  be  good,  one  must  be  rather  cautious  in 
calibrating  a  detector  of  singlet  oxygen  solely  on  its  reported  performance 
and  efficiency,  as  has  been  sxiggested  (Ref.  4). 

The  emission  resulting  from  NO  0  recombination  is  broad  band  and  partially 
overlaps  the  narrow  O^(a-X)  emission  at  1.25  to  1.29  fim.  By  using  two  narrow 
band-pass  filters  at  1.268  ±  0.010  pm  and  1.3152  ±  0.010  pm.  the  ratio  of  the 
observed  intensities  for  NO  0  and  O^ta-X)  emission  may  be  experimentally 
determined.  The  NO  -t-  0  emission  at  1.3152  pm  was  undetectable  under  typical 
operating  conditions  and  a  massive  air  leak  was  introduced  for  these  tests. 

For  an  air  flow  rate  upstream  of  the  microwave  generator  of  m  600  seem  and  a 
total  pressure  of  1  torr,  a  green  emission  was  visually  observed  in  a  darkened 
room  and  a  0.3  mV  signal  recorded  at  1.315  pm  on  the  IR  detector.  The  artifl- 

3  4 

cial  leak  rate  is  10  to  10  times  larger  than  the  apparatus  leak  rate.  Most 
of  the  apparatus  leak  rate  occurs  downstream  of  the  microwave  discharge. 

Thus,  an  upper  bound  to  the  NO  -«■  0  signal  under  typical  operating  conditions 
is  0.3  pV,  which  is  100  times  smaller  than  the  detector  dark  signal. 

In  the  absence  of  a  mercury  coating  on  the  exit  walls  of  the  discharge  tube, 
oxygen  atoms  were  detected  at  the  EPR  cavity.  With  small  amounts  of  mercury 
(Hg)  applied  to  the  walls,  the  oxygen  atom  concentration  was  reduced  by  50 
percent,  and  with  moderate  mercury  coating  no  oxygen  atoms  were  detected  at 
the  EPR  cavity. 


it 
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3.3.2  Optical  Collection  and  Detection 

The  IR  emission  detection  system  including  detector  and  electronics  exhibits 
a  highly  linear  response  to  incident  Intensity,  as  shown  in  Figure  18.  The 
emission  from  a  blackbody  source  was  transmitted  to  the  collimating  lens, 
fiber  optic,  and  germanium  detector  through  an  interchangeable  neutral  density 
filter  and  the  resultant  AC  signal  was  recorded  on  the  lock-in  amplifier.  The 
observed  signals  range  from  0.2  to  300  mV  and  are  nearly  linear  over  the  full 
range  with  a  correlation  of  >  0.9994.  Thus,  the  EPR  calibration  can  be 
extrapolated  to  considerably  larger  O^Ca)  concentrations  with  high  confidence. 
Note  that  with  the  calibration  factor  presented  in  Section  4.0,  the  300  mV 
signal  corresponds  to  an  0^(a)  concentration  of  12.7  torr. 

The  data  presented  in  Figure  18  determine  the  relative  calibration  of  the  two 
different  IR  detectors.  The  relative  slopes  are  eqxial  to  the  relative  cali¬ 
brations  as  follows: 

CALj/CAL^  ■  slope^/slope^j  (27) 

where: 

slope^  ■  slope  for  detector  i 

CALj  m  calibration  factor  for  detector  i,  as  defined  in  Eqxiatlon  27 

The  ratio  CAL^/CAL^  is  determined  to  within  2  percent  amd  is  considerably  more 
accurate  than  any  absolute  calibration.  This  result  can  be  used  to  further 
constrain  the  error  of  two  independent,  absolute  calibrations  and  to  transfer 
a  known  calibration  from  one  detector  to  a  second. 

The  IR  detector  response  is  strongly  dependent  on  chopping  frequency.  The 
teqptoral  response  of  the  ADC  model  403L  germanium  detector  is  relatively  slow 
(10  ms).  The  ADC  detector  responds  to  light  modulated  in  the  range  of  0.5  to 
500  Hz.  The  frequency  stability  of  the  mechanical  light  choppers  was  within 
±0.2  percent. 
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Characterizing  the  collination  length  of  the  optical  collection  system  is  a 
difficult  task.  Several  qualitative  checks  were  performed  to  assess  the 
sensitivity  of  the  calibration  to  viewing  volume.  First,  the  collection  optic 
was  moved  away  from  the  fixed  calibration  duct  along  the  optical  axis  by  0  to 
2  cm  and  no  change  in  signal  was  observed.  Second,  the  O^(a-X)  emission  was 
observed  at  both  detection  blocks  B3  and  B2  (path  length  2.21  cm)  and  at  the 
calibration  duct  (path  length  3.99  cm).  Comparisons  between  the  signals  at 
B3,  B2  and  at  the  duct  were  made.  Corrections  for  the  path  lengths  and  the 
pressure  drops  were  made,  and  an  Independent  estimate  of  the  deactivation  was 
provided  by  DEACT  code  (Section  4.0).  The  values  measured  for  02(a-X)  at  B3 
and  the  corrected  value  at  the  duct  agreed  to  within  2  percent.  Similar 
results  were  obtained  for  the  values  at  the  duct  amd  B2.  Finally,  a  large 
path  length  correction  of  4.8  was  successfully  applied  to  actual  COIL  hardware 
between  the  oxygen  generator  and  the  diagnostic  duct.  Thus,  viewing  volume 
corrections  based  on  relative  optical  path  lengths  is  viable  for  volume  ratios 
of  2  to  5.  It  is  still  highly  desirable  to  calibrate  on  an  exact  volume 
replica  of  the  intended  laser  hardware  for  the  best  calibrations. 

3.3.3  The  EPR  Spectrometer 

The  EPR  spectrometer’s  signal  to  noise  characteristic  was  measured  with  a 
standard  pitch  sample  (Ref.  40).  The  noise  trace  was  obtained  with  no  pitch 
sample  in  the  cavity.  The  signal  to  noise  ratio  obtained  from  the  measured 
spectrum  is  445:1,  in  good  agreement  with  the  instrumental  specification  of 
500:1. 

The  0_(X^  ")  K  =  3,  J  *  4,  M  «  3  4  transition  at  H  ®  9188  G  was  used  for 

the  calibration  reference  in  the  present  study.  The  spectriim  is  dense  below 
a  9000  G  and  no  attempt  to  assign  these  transitions  was  made. 

A  complete  EPR  spectrum  of  O^(a^A)  J  »  2  is  shown  in  Figure  19.  The  four 
AM  >  1  transitions  are  indicated  and  are  separated  by  97  to  94  G.  The 
relative  Intensities  are  in  the  proportion  2.09  :  3.11  :  2.99  :  1.91,  in  good 
agreement  with  theory.  The  strong  feature  at  H  ■  9673  G  is  due  to  the  02(X^) 


41 


bitegrated  E8R  signal  Obssrved  E8R  signal 

(arbitrary  units)  «  (arWtrary  units) 


»70  9570  9770 
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a)  Observed  ESR  signal  of  0^(a) . 


9370  9570  9770 


H  (gauss) 

b)  First  integral  of  observed  signal  showing  absorption  feature. 
Figure  19.  Complete  ESR  spectrum  of  02(a)  J  ■  2  at  v  -  8.92  GHz. 
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transition.  The  0  (a^A)  M  »  -1  ^  0  and  -2  ->  -1  transitions  should  be  avoided 

2  j 

In  these  calibration  procedvres  due  to  convolution  effects  with  this  strong 
ground  state  feature  (Fig.  19).  Thus,  the  stronger  of  the  remaining  lines, 
the  M  a>  0  1  transition  at  H  »  9522  G,  was  typically  used  to  monitor  0  (a^A) 

J  2 

concentrations.  A  high  resolution  scan  of  the  0^(a  A)Jb2M  «0->1  transi- 

2  J 

tion  centered  at  H  ®  9522  G  is  shown  in  Figure  20  and  the  O^(X^)  K  »  3  J  ■  4 
«  3  4  transition  at  H  =  9188  G  is  shown  in  Figure  6.  A  strong  signal 

with  moderate  modulation  is  achieved  and  a  nearly  complete  integration  of  the 
signal  over  the  20  G  scan  width  is  obtained. 

The  tramsitions  shown  in  Figure  19  are  pressure  and  modulation  broadened. 
Characterization  of  the  line  shape  profiles  including  broadening  coefficients, 
Voight  fits,  and  spectral  area  integrations  is  described  in  the  following 
paragraphs. 

Pressure  broadening  has  a  significant  effect  on  the  line  shapes  observed 
during  the  calibration  procedure.  The  pressure  broadening  of  the  0  (X^>  by 
collisions  with  oxygen  has  been  previously  reported  as  1.40  (Ref.  43)  to  1.87 
G/Torr.*  The  pressure  broadening  observed  in  the  present  apparatus  is  shown 
in  Figure  21  and  provides  a  broadening  coefficient  of  1.70  ±  0.16  G/Torr.  The 
observed  line  shapes  are  also  modulation  broadened.  It  has  been  previously 
established  both  experimentally  and  theoretically  that  the  area,  or  second 
integral,  of  an  EPR  spectrum  increases  linearly  with  modulation  amplitude 
(Ref.  43).  As  seen  in  Figure  22a  this  linear  relationship  is  verified  \inder 
the  operating  conditions  of  this  work.  The  peak-to-peak  width  of  the  EPR 
signal  also  increases  with  modulation  amplitude  (Fig.  22b). 

Modulation  broadening  has  been  examined  in  detail  (Refs.  28-32).  Modulation 
broadening  generally  affects  the  central  portion  of  the  EPR  spectrum  and  the 
wings  are  relatively  unaffected  (Ref.  23).  For  a  Lorentzian  profile  with 
modulation  broadening,  the  maximum  EPR  signal  occurs  idien  the  modulation 


*  Ibid,  page  17. 
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Magnetic  field  atrength  (gauaa) 


a)  Observed  ESR  signal. 


b)  First  Integral  of  ESR  signal  (absorption  profile) . 


c)  Second  Integral  of  observed  signal  used  for  obtaining  areas. 


Figure  20.  Hl^  resolution  ESR  spectrum  of  0  (a)  J  -  2,  M  -  0->l  transition 

2  *7 

at  >/  -  8.92  GHz. 
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Figure  21,  Pressure  dependence  of  observed  0^(X)  ESR  line  shape. 
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anplitude  is  be  1.5  times  the  peak-to-peak  width  (Ref.  23).  Modulations  this 
great  could  not  be  reached  with  the  apparatus  used  in  this  work,  but  the 
observed  areas  increase  to  modulations  as  high  as  75  percent  of  the  peak-to- 
peak  width.  Figure  23  shows  the  observed  line  shapes  are  strongly  affected 
for  modulations  >  50  percent  of  the  line  width.  Typically,  a  modulation 
intensity  of  H  =  0.8  G  was  used  to  avoid  these  strong  modulation  effects. 

The  observed  line  shapes  were  fit  to  Lorentzian,  Gaussian,  and  Voight  profiles 
using  a  Gauss-Newton  routine.  Figure  24  shows  the  comparison  of  these  fits 
with  the  observed  data.  The  observed  line  shape  is  intermediate  between 
(kiussian  and  Lorentzian  and  is  reasonably  well  fit  by  the  convolved  Voight 
profile.  The  Voight  profile  for  0^(x\)  underestimates  the  signal  far  from 
resonance  by  as  much  as  20  percent,  but  reproduces  the  integrated  area  to 
within  5  percent. 

Similar  results  are  achieved  for  O^(a^A)  with  the  Voight  fit  shown  in  Figure 
24b.  For  an  at  a  total  pressure  of  0.48  torr,  at  the  EPR  the  ratio  of 

Gaussian  width  to  Lorentzian  width,  t)  «  Av/(Av/in2),  is  i)  >=  3.8,  indicating 

d  L 

the  Lorentzian  contribution  to  the  line  width  is  m  27  percent  of  the  resultant 
width  and  the  Gaussian  contribution  is  m  86  percent  (Ref.  27). 

The  relative  intensities  of  the  02(a)  and  ©^(X)  EPR  transitions  were  required 
for  the  calibration  (Eq.  8).  An  observed  EPR  spectrum  was  proportional  to  the 
derivative  of  the  absorption  profile,  thus  a  double  integration  of  the  spectra 
led  to  the  relative  intensities.  The  integration  was  limited  to  a  finite  mag¬ 
netic  field  scan  due  to  interferences  from  other  spectral  features.  In  the 
present  study  scan  widths  of  20  or  40  G  were  typically  used  and  constituted  5 
to  10  times  the  width  of  the  spectral  feature.  Figure  25  illustrates  the 
error  associated  with  truncating  the  integration  for  both  Gaussian  and  Lorent¬ 
zian  line  shapes.  The  scan  widths  of  the  present  study  introduces  essentially 
no  error  for  Gaussian  profiles,  but  as  much  as  30  percent  error  for  Lorentzian 
profiles  (Ref.  23). 

Since  the  observed  profiles  are  only  approximately  fit  by  a  Voight  line  shape, 
an  experimental  test  of  the  truncation  error  was  required.  A  Voight  fit  to  a 
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(H«/AHp_p)/2 


Figure  23.  Noinnallzed  ESR  profile  width  as  a  function  of  magnetic  field 
modulation  amplitude. 
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a)  Lorentzian  and  gaussian  fits  to  the  integrated  ESR  profile. 


b)  Voight  fit  to  integrated  ESR  profile. 


Figure  24.  Comparison  of  integrated  ESR  profiles  with  standard  line  shapes. 
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20  G  scan  of  the  0  (X)  line  used  the  calibration  was  made.  Using  these  fit 

2 

parameters,  the  Volght  profile  was  extended  to  a  width  of  up  to  200  G.  Figure 
26  presents  the  dependence  of  the  calculated  integrated  area  for  various  field 
scan  widths  of  these  Voight  profiles.  The  integrated  area  versus  scan  width 
is  relatively  flat,  but  does  increase  by  about  6  percent  from  a  20  G  scan  to  a 
200  G  scan.  This  error  is  intermediate  between  the  Lorentzian  and  Gaussian 
errors  shown  in  Figure  25.  The  double  integral  of  the  observed  absorption 
feature  also  supports  this  near  complete  integration.  Figure  20c  shows  the 
spectral  area  as  a  function  of  field  strength  and  is  nearly  flat  for  field 
widths  above  25  to  30  G. 

No  theoretical  correction  for  the  integral  truncation  error  was  made  in  the 
present  study.  A  Lorentzian  wing  correction  was  not  Justified  since  the  line 
shape  had  considerable  Gaussian  nature.  The  truncation  error  was  small  for  a 
given  spectrum  and  was  further  reduced  since  only  the  ratio  of  two  observed 
areas  was  required  for  the  detector  calibration.  A  quantitative  evaluation  of 
the  truncation  error  is  provided  in  Section  5.0. 

An  inhomogeneous  magnetic  field  may  lead  to  asymmetric  ESR  absorption  spectra. 
The  symmetry  of  an  ESR  profile  may  be  determined  from  the  first  integral  of 
the  original  ESR  profile.  A  symmetric  absorption  profile  has  even  symmetry 
about  the  profile  peak.  The  symmetry  of  the  first  integral  profile  was 
checked  by  comparing  pairs  of  points  equidistant  in  field  from  the  peak.  In 
the  present  experiment  the  maximum  difference  between  a  pair  of  points  was  1.2 
percent,  indicating  near  symmetry.  This  symmetry  Justified  averaging  the 
points  in  an  observed  ESR  spectrum  to  calculate  the  baseline  (Ref.  23). 

The  effects  of  microwave  power  saturation  were  examined  in  the  present  study 
and  the  results  are  presented  in  Figure  27.  If  the  rate  of  stimulating 
transitions  between  spin  states  approaches  the  rate  of  spin  state  relaucation, 
the  transition  becomes  power  saturated  (Ref.  23).  Saturation  has  been 
observed  on  the  02(x\)  K  =  3,  J  *  4,  »  3  -»  4  transition  at  microwave 

powers  of  about  40  mU  (Ref.  6).  A  linear  fit  to  the  low  power,  <100  mW, 
portion  of  the  curve  provides  a  near-zero  intercept  and  a  good  correlation  of 
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r  B  0.973.  The  curve  of  Figure  27  begins  to  deviate  froa  this  linear  response 
at  Klystron  powers  above  100  mW.  To  avoid  the  possible  cosplicatlons  of  satu¬ 
ration  effects,  microwave  powers  of  40  mU  were  used  in  the  present  studies. 

The  relative  EPR  gain  settings,  G^.  are  critical  to  the  absolute  O^Ca^A) 
concentration  calibrations,  as  seen  in  Eq\iation  10.  The  relative  peak 
intensities  are  plotted  as  a  function  of  the  gain  up  to  8000  in  Figure  28. 

The  y-intercept  is  zero  within  error,  -0.15  ±  0.19,  and  the  correlation  to  a 
linear  fit  is  0.99987.  The  Instrumental  gain  settings  were  accurate. 

The  linearity  of  the  EPR  apparatus  was  verified  by  observing  the  integrated 
areas  of  the  0^(X)  EPR  absorption  features  for  various  flow  rates  of  pure 
0^(X).  Figure  29  demonstrates  a  linear  EPR  response  with  a  correlation  of 
r  B  0.9933  and  a  standard  deviation  of  O'  ^  30.8.  The  intercept  is  zero 

within  the  uncertainty  of  the  data.  The  scatter  in  the  data  of  Flgvire  29  is 
primarily  due  to  the  reproducibility  of  observed  EPR  integrated  absorption 
areas . 

Typically,  the  ratio  of  areas  from  0  (^A)  and  0  (X^)  spectra  were  calculated 
to  determine  the  fraction  of  A)  to  total  oxygen  concentration.  This 
assumed  that  with  the  microwave  discharge  on,  the  oxygen  flow  consisted  of 
only  these  two  species.  Another  method  for  calculating  the  ratio  of  O^('a)  to 
total  oxygen  which  makes  this  assumption  is  to  compare  the  EPR  signals  of 
©^(X^r)  with  the  microwave  discharge  on  with  the  microwave  off. 

l0,(a*A)l  =  [0,(X^Z)1  -  [0,(X^Z)]  (28) 

2  2  o  2 

where: 

I0^(X  Z)]  ■  ground  state  concentration  when  microwave  is  off 

2  o 

[O^(XZ)]  ®  ground  state  concentration  when  microwave  is  on 

This  difference  method  is  not  preferred  since  small  changes  in  large  signals 

are  obtained.  Figure  30  Illustrates  the  ratio  (0  (^A)l/[0  (X^)l,  calculated 

2  2 

by  the  difference  method  (A)  and  the  observed  0^{a)  EPR  signals  (+).  The  two 
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Figure  29.  Area  under  absorption  profile  as  a  function  of  0^(X)  pressure. 
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Figure  30.  Comparison  between  difference  and  ratio  methods  for  determining 
the  fraction  of  0  (a) . 
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■ethods  agree  to  within  2  percent.  This  approach  provided  a  quick  system 
check  and  indicated  that  the  microwave  discharge  produced  no  high  concentra¬ 
tions  of  other  excited  states  of  oxygen. 

The  difference  in  these  two  results  could  be  attributed  to  the  presence  of 

0  (b^£).  An  estimate  for  the  concentration  of  0  (b)  can  be  made  from  the 
2  2 

kinetic  analysis  presented  in  Section  2.4.  The  steady-state  concentration  of 
O^Cb)  is 

tO,(^S)]  =  k  ,[0,(^A)]V(k^  +  k  ^tO,(^Z)])  (29) 

2  13  2  20  18  2 

For  a  wall  rate  of  k  -20  s~^  which  corresponds  to  a  deactivation  probability 

20 

-2  3  -1 

of  y  »  10  in  a  2.S4-cm  diam  tube  at  a  velocity  of  10  cm  s  ,  the  singlet 

12 

sigma  concentration  is  a  10  molecules  per  cubic  centimeter,  or  0.1  percent 

of  the  the  O^(^A)  concentration.  This  calculation  is  sensitive  to  the  wall 
2 

rate  and  if  a  slower  rate  is  chosen  so  that  the  deactivation  is  controlled  by 
gas  phase  quenching,  then  an  upper  bound  to  the  0^(b)  concentration  of  ~  5 
percent  of  [O^Ca^A)]  is  obtained.  Thus,  the  differences  obser\  a  in  Figure  30 
could  be  partially  attributed  to  0^(b)  or  may  simply  reflect  systematic 
differences  in  the  measurements. 
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4.0  RESULTS  AND  CALIBRATIC»4  DATA 


The  infrared  detection  system  calibration  is  determined  by  plotting  the 
concentration  of  0  (a^A)  versus  the  observed  detector  voltage,  V  . 

2  IR 

mCa^A)]  =  CAL  (30) 

The  detector  voltage  is  measured  directly  at  the  diagnostic  duct.  The  O^(a^A) 
concentration  is  calculated  from  the  observed  EPR  spectra  and  corrected  for 
transport  losses  to  provide  the  concentration  at  the  diagnostic  duct.  The 
analysis  of  Section  2.4  is  used  below  to  derive  the  quantity  [0^(a^A)l  at  the 
diagnostic  duct. 


Equation  31  provides  the  relationship  between  cbserved  EPR  spectra  and 
concentration  of  O^ta^A)  at  the  EPR  cavity: 

lO,(a‘A)l^„„  »  [r/(l  +  r)l  (P„„/K  T)  (31) 

2  EPR  EPR  B 


where: 


'■o  ,.,''■0  ,x,>  «=o  ,.,^=0  <x,> 
2  2  2  2 


The  pressure  at  the  EPR  cavity,  P  .  is  calculated  from  the  friction  loss 
Equation  33  and  the  observed  pressures  at  detection  blocks  B1  and  B2.  The 
distances  are  measured  from  the  EPR. 


P  ^  «  P  ^  + 

EPR  1 


(P,^  -  p/)  (X  -  X  )  /  (X  -  X  ) 
2  1  EPR  1  2  1 


where: 

x^  a  position  of  detection  block  B2  =  41.9  cm 

x^  ■  position  of  detection  block  B1  =  54.6  cm 

X  ■  position  of  EPR  cavity  =  0  cm 
P^m  pressure  at  block  B1 

Pg"  pressure  at  block  B2 


(33) 
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The  gas  temperature  In  Equation  8  was  measured  at  blocks  B3,  B2,  and  at  the 
duct  and  was  Independent  of  position,  x.  The  temperature  at  the  calibration 
duct  was  measured  8.9  cm  from  the  optical  windows. 

The  data  from  a  typical  calibration  conducted  on  19  February  1987  are  provided 
in  Table  4.  For  the  locations  of  various  measurements,  refer  to  Figure  10. 
Reducing  the  data  according  to  the  above  procedure  provides  the  calibration 
plot  shown  in  Figure  31.  The  plot  had  a  linear  correlation  of  r  -  0.9980,  a 
near  zero  intercept  of  0.002  ±  0.01  and  a  slope  calibration  factor,  CAL  »  1.60 
X  lo'®  ±  0.01  cm'^. 


Table  4.  Sample  calibration  data 

Run 

No.  Pressure  (torr)  Signal  (volts)  Flow  rate  (seem) 

B3  Duct 


B1 

B2 

B3  Duct 

B1 

1 

0.402 

0.526 

0.575 

0.554 

0.55 

2 

0.264 

0.346 

0.380 

0.365 

0.33 

3 

0.163 

0.214 

0.236 

0.226 

0.23 

4 

0.098 

0.127 

0.141 

0.135 

0.11 

5 

0.345 

0.453 

0.495 

0.477 

0.44 

6 

0.405 

0.531 

0.582 

0.560 

0.52 

7 

0.511 

0.669 

0.731 

0.705 

0.60 

8 

0.593 

0.775 

0.847 

0.817 

0.67 

9 

0.742 

0.966 

1.055 

1.02 

0.73 

10 

0.862 

1.12 

1.22 

1.18 

0.79 

(Cable) 

A 

B 

C 

0.72 

1.3 

1.0 

1.4 

1.22 

201 

0.46 

0.91 

0.75 

1.07 

0.95 

201 

0.32 

0.75 

0.46 

0.65 

0.57 

52 

0.16 

0.52 

0.24 

0.34 

0.32 

30 

0.62 

1.13 

0.85 

1.16 

1.04 

149 

0.70 

1.25 

0.95 

1.29 

1.18 

200 

0.82 

1.43 

1.17 

1.60 

1.38 

300 

0.90 

1.53 

1.24 

1.73 

1.52 

402 

1.0 

1.67 

1.36 

1.92 

1.65 

602 

1.08 

1.75 

1.46 

2.08 

1.80 

800 

This  calibration  does  not  account  for  02(a*A)  deactivation  during  the  trans¬ 
port  from  the  diagnostic  duct  to  the  EPR  cavity.  The  analysis  of  Section  2.4 
provides  the  basis  for  making  the  necessary  deactivation  correction.  The 
transport  of  O^(a'A)  is  described  by  Equation  34: 

[0,(a*A)]“^  =  {  (0,(a^A)]  "^  +  A/B  }  exp  (Bt)  -  A/B  (34) 

2  2  o 
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The  deactivation  observed  in  the  present  apparatus  is  snail,  usually  less  than 
20  percent  between  block  B3  euid  the  EPR.  and  Equation  34  can  be  simplified  by 
the  approximation  exp(Bt)  a  1  -t-  Bt: 

[0  (a^A)!'^  *  (0,(a^A)]  {  m  (a^A)]  '^  +  Ay^  >  Bt  (35) 

2  2  o  2  o 

This  inverse  linear  relationship  provides  for  a  simple  deactivation  correc¬ 
tion.  In  Figure  32.  the  concentration  of  0^(a)  is  shown  as  a  fiuiction  of 
flow  distemce. 

The  IR  detector  signals  at  blocks  Bl,  B2,  and  B3  directly  provide  a  relative 

measure  of  the  0  (a^A)  concentration  and  were  used  to  assess  the  value  of  the 
2 

term: 


(  [O(a^A)]"^  +  }  Bt  (36) 

That  is,  the  slope  is  determined  by  the  measured  IR  detector  signals  and  a 
linear  Interpolation  is  used  to  correct  for  deactivation  between  the  diagnos¬ 
tic  duct  and  the  EPR  cavity.  Note  that  the  observed  IR  detector  signal 
depends  on  both  the  fraction  of  O^(a^A)  in  the  flow  and  the  total  pressure. 

The  pressure  drops  between  detection  blocks  due  to  friction  losses  and  must  be 
accounted  for  by  the  relationship: 

«  P  (  [0,(a^A)l  /  [N^]  )  (37) 

m  2  X 

Thus,  the  correct  dependent  variable  for  Figure  31  is: 

[0,(a‘A)]‘'  “  P  /  C38) 

Since  the  diagnostic  duct  and  entrance  ramps  are  spatially  symmetric  and  the 
optical  viewing  port  is  at  the  midpoint,  one-half  of  the  observed  deactivation 
between  blocks  B3  and  B2  is  attributed  to  the  deactivation  between  the  diag¬ 
nostic  duct  and  block  B2.  Almost  no  deactivation  occurs  between  blocks  B2 
and  Bl  due  to  high  flow  velocity  in  this  narrow  channel. 
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The  calibration  plot  corrected  for  deactivation  by  the  above  procedure  Is 
given  in  Figure  31.  Again,  the  plot  Is  quite  linear  with  a  correlation  of 
r  “  0.9972,  has  a  near  zero  Intercept,  -0.05  ±  0.01,  and  accurately  defines 

15  —3 

the  calibration  factor,  CAL  =  1.41  ±  0.01  x  10  cm  The  transport  deacti¬ 
vation  correction  reduces  the  calibration  by  a  5  percent. 

The  calibration  must  provide  a  linear  plot  with  zero  Intercept.  Otherwise, 
significant  systematic  errors  are  Indicated.  In  particular,  a  nonzero  Inter¬ 
cept  may  result  from  significant  atmospheric  leaks  and  a  high  background 
pressure . 

A  more  sophisticated  study  of  the  kinetic  deactivation  due  to  transport  was 
conducted  using  the  DEACT  computer  model  (Ref.  42).  A  complete  description  of 
the  gas  phase  kinetics  of  oxygen  CRef.  33)  was  coupled  with  a  two-dimensional 
fluid  dynamics  code  to  predict  O^Ca^A)  concentrations  as  a  function  of  flow 
distance  for  the  calibration  apparatus  geometry.  A  comparison  of  the  observed 
data  and  model  predictions  is  shown  in  Figure  32.  Deactivation  is  greatest  in 
the  diagnostic  duct  where  the  flow  velocity  is  lowest.  Also  shown  in  Figure 
32  is  a  linear  interpolation  between  blocks  B3  and  B2.  While  the  deactivation 
Is  not  linear  bet;;3en  these  blocks.  It  Is  symmetric  and  the  linear  and  exact 
calculations  agree  at  the  midpoint. 

The  data  of  Figure  32  were  not  well  fit  by  a  single  deactivation  rate.  The 
modeled  results  were  obtained  with  a  single  parameter,  the  wall  deactivation 
rate,  defined  as  a  probability,  y.  To  successfully  reproduce  the  data,  the 
parameter  y  was  allowed  to  vary  as  a  function  of  flow  rate.  While  this 
approach  is  not  physical,  it  does  allow  for  good  agreement  with  the  data.  A 
more  precise  and  detailed  study  to  characterize  the  deactivation  of  singlet 
oxygen  on  surfaces  is  underway  (Ref.  44). 
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5.0  E3W0R  ANALYSIS 

The  purpose  of  this  section  Is  to  address  all  known  sources  of  error  associa¬ 
ted  with  the  calibration.  No  rigorous  statistical  analysis  is  presented; 
however,  the  estimated  error  in  measured  quantities  is  given.  Considerable 
thought  has  gone  into  addressing  all  possible  sources  of  systematic  errors. 
The  most  important  conclusion  of  this  section  is  that  the  single  largest 
source  of  error  in  the  calibration  is  the  uncertainty  associated  with  the 
singlet  delta  EPR  signal.  Refer  to  List  Of  Symbols  for  definitions  of  the 
variables  shown  in  the  tables. 

5.1  STATISTICAL  ERRORS  OF  MEASURED  AND  DERIVED  QUANTITIES 

The  random  (or  statistical)  error  in  each  measurable  quantity  has  been 
estimated,  based  on  experience  gained  from  performing  hundreds  of  EPR  scans. 
These  estimated  errors  are  summarized  in  Table  5.  Since  the  measured  values 
nearly  always  fell  within  these  error  limits,  these  errors  can  be  considered 
estimates  of  the  95  percent  confidence  limit. 


Table  5. 

Resolution  and 

statistical  error 

in  observed  quau\tities. 

Measurable 

Quantity 

Instrument 

Resolution 

Typical  Value 

Estimated  95  Percent 
Confidence  Limit 

T 

0.1  K 

294  K 

0.3% 

P 

0.2  % 

0.1  -  1.5  Torr 

0.5% 

m 

1.0  •/. 

50  -  2000  seem 

0.5% 

Ax 

1  mm 

1000  mm 

0.1% 

r 

o 

0.1  mm 

10  mm 

1 . 0% 

V 

IR 

0.01  mV 

0.3  -  2.0  mV 

1 . 0% 
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The  errors  in  derived  quantities  such  as  velocity,  lO^(a'A)],  and  the  calibra¬ 
tion  factor,  Cal,  are  obtained  by  the  appropriate  propagation  of  errors  (Ref. 
45)  and  are  reported  in  Table  6.  Briefly,  the  error  in  a  derived  quantity  is 
the  total  derivative  of  that  quantity  divided  by  the  value  of  the  quantity, 
dq/q,  where  the  total  derivative  is  obtained  by  partial  differentiation  with 
respect  to  all  the  independent  variables.  In  some  cases  the  analytic  expres¬ 
sion  for  the  partial  derivative  has  been  approximated  to  provide  a  constant, 
numerical  coefficient  to  the  error  in  measured  quauritities.  The  derived  errors 
are  provided  in  Table  6. 

By  far  the  greatest  source  of  statistical  error  is  provided  by  the  observed 
EPR  spectrum,  a'(H),  for  O^Ca).  Indeed,  the  12  percent  variations  in  observed 
EPR  absorption  signals,  due  to  the  poor  signal-to-noise  ratios  at  these  low 
0^(a)  concentrations,  directly  translates  into  a  20  percent  error  in  the 
resultant  calibration.  None  of  the  observed  errors  are  limited  by  Instrument 
resolution  or  the  precision  of  required  constants. 


Table  6.  Statistical  errors  of  derived  quantities. 


Quantity 

Calculated 

Error 

pressure  in  EPR  cavity  (P  ) 

EPR 

0.9% 

gas  velocity  (v) 

3.8% 

Integral  of  Absorption  (I  ) 

1 . 0% 

Integral  of  Absorption  (I^ 

2 

12.0% 

Fraction  of  singlet  delta  by  ratio 

method 

13.0% 

Fraction  of  singlet  delta  by  difference  method  18.0% 

Partial  pressure  of  Singlet  Delta 

2 

14.0% 
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5.2  SOURCES  OF  SYSTEMATIC  ETOORS 

The  potential  sources  of  systematic  error  are  identified  in  Table  7.  Extreme 
care  was  taken  to  check  for  systematic  errors  and  eliminate  them.  Following 
the  table  is  a  discussion  of  each  source  of  systematic  error. 

Table  7.  Primary  sources/symptoms  of  systematic  errors. 

nonzero  Cal  plot  intercept 

nonlinear  Cal  plot 

nonzero  pump  curve  intercept 

Assumption  that  only  O^CX.a)  present 

integral  truncation 

spectral  interference 

nonsteady  flow 

deactivation  correction 

improper  H  amplitude  calibration 

EPR  signal  (a'  (H)) 

EPR  signal  (a'  (H) ) 

0  a 
2 

EPR  Gain  Setting  (G) 

Magnetic  Field  Scan  (AH  ) 

scan 

Modulation  (H  ) 

n 

EPR  Microwave  Power  (P^) 


A  nonzero  intercept  in  a  calibration  plot  is  a  strong  indication  of  serious 
systematic  errors.  Experimental  problems  such  as:  (1)  the  detection  of 
background  light  or  extraneous  emissions,  (2)  substantial  vacuum  leaks,  (3) 
significant  kinetic  deactivation,  (4)  nonlinear  response  of  detection  systems 
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(IR  sensor  or  EPR  spectrometer),  or  (5)  equipment  malfunction  can  lead  to  such 
non-zero  intercepts.  As  shown  in  Figure  31,  the  intercepts  observed  in  the 
present  study  are  zero,  within  experimental  error.  This  should  be  the  case 
for  any  reliable  calibration.  Similarly  the  calibration  plot  should  exhibit 
no  curvature.  The  linear  correlation  coefficient  for  the  calibration  plot 
corrected  for  deactivation  is  r  =  0.9979,  and  thus  no  significant  systematic 
errors  were  apparent. 

A  large  leak  rate  or  zero  offset  in  mass  flow  controller  or  pressure  trans¬ 
ducer  can  lead  to  a  nonzero  intercept  in  Figure  18.  While  the  calibration 
factor  does  not  depend  on  mass  flow  rate  directly,  the  deactivation  correction 
does  depend  on  velocity  which  is  obtained  from  flow  rates.  This  constant 
offset  effect  will  cause  the  greatest  percentage  error  at  very  low  flow  rates 
and  longest  deactivation  times.  Since  the  deactivation  correction  is  always 
less  than  12  percent  in  the  present  apparatus  and  an  offset  of  5  percent  of 
the  total  flow  will  cause  a  small  error  in  velocity,  and  a  negligible  error  in 
the  calibration. 

The  presence  of  species  other  than  0^(X)  and  0^(a)  in  the  flow  from  the 
microwave  source  is  rather  limited.  Typically  the  concentration  of  O^(b'z)  is 
about  1  percent  of  the  total  0^  flow  and  oxygen  atoms  account  for  only  about  1 
in  10  particles.  These  results  agree  with  the  predictions  of  the  difference 
method  for  determining  the  concentration  of  0^(a)  provided  by  Equation  8.  The 
error  in  integrated  area  introduced  by  the  finite  field  scan  width  of  20  G  was 
established  in  Section  3.0  as  <  6  percent. 

The  truncation  error  always  underestimates  the  area  of  the  spectral  feature. 
Since  only  the  ratio  of  0  (a)  concentration  with  respect  to  0  (X)  is  needed, 
it  is  the  difference  in  truncation  errors  for  the  two  scans  that  introduces 
any  systematic  error.  It  is  more  difficult  to  experimentally  determine  the 
truncation  error  for  02(a)  due  to  the  greater  extent  of  noise  in  the  spectrum. 
Thus  the  relative  truncation  error  can  only  be  estimated.  The  truncation 
error  is  primarily  dependent  on  the  ratio  of  spectral  feature  width  (AH  )  to 

pp 

the  scan  width  (AH  ).  The  absorption  profile  widths  of  0  (a)  and  0  (X) 

•can  22 


65 


WL-TR-90-85 


were  s  the  same  within  experimental  error.  Therefore,  the  relative  truncation 
error  is  negligible. 

No  emission  signals  were  observed  at  the  IR  detector  above  the  dark  signal  of 
0.01  mV.  Neither  background  radiation  from  the  room  nor  NO  +  0  recombination 
emission  were  detectable  in  the  present  apparatus  and  introduced  no  error  in 
the  calibration. 

The  microwave  source  of  0^(a)  must  be  extremely  stable  to  prevent  any  changes 
in  flow  conditions  or  partial  pressures  during  the  time  required  to  obtain  all 
EPR  spectra.  The  characterization  of  the  present  apparatus  established  no 
long  term  drift  in  the  0^(a]  source  and  local  fluctuations  of  <0.5  percent.  A 
chemical  source  of  0^(a)  would  not  be  nearly  as  stable. 

The  observed  EPR  line  positions  and  line  widths  agree  well  with  those 
previously  reported.  However,  no  direct  determination  of  the  magnetic  field 
strength  was  accomplished  in  the  present  study.  The  present  calibration 
facility  should  be  modified  to  provide  a  direct  determination  of  magnetic 
field  strength. 

Proper  maintenance  of  the  fiber  optics  bundles  is  critical.  The  detector  col¬ 
lection  efficiency  can  be  severely  affected  by:  (1)  changing  or  re-installing 
the  fiber  optic  bundle,  (2)  stressing  or  bending  the  fiber  bundle,  (3) 
changing  the  very  sensitive  focus  of  the  collimating  optics  by  flexing  the 
cable  or  re-orienting  the  optic  mount,  or  (4)  chipped  or  dirty  fiber  end 
surfaces.  The  complete  detection  system  including  collection  optics  and  fiber 
optics  should  be  maintained  as  an  integral  system. 

The  infrared  detectors  must  be  maintained  at  liquid  nitrogen  temperatures  to 
insure  maximum  sensitivity  since  the  detector  response  is  temperature- 
sensitive.  Also,  the  voltage  supplied  to  the  detector  must  remain  at  -  250  V 
DC  for  proper  operation. 
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In  applying  the  diagnostic  to  COIL  devices,  one  must  be  careful  to  ensure  that 
the  windows  or  diagnostic  ports  on  the  device  are  not  affected  by  device  oper- 
tions.  Often,  the  windows  can  be  coated  with  salts  or  other  flow  contaminates 
that  may  severely  degrade  the  detector  collection  efficiency. 

Periodic  checks,  possibly  through  relative  blackbody  calibrations,  are  criti¬ 
cal  to  assessing  any  changes  in  the  performance  of  the  diagnostic  during  a 
series  of  laser  test  sequences. 

5.3  TOTAL  ERRORS  ASSOCIATED  WITH  CALIBRATION 


In  light  of  the  previous  discussions,  the  following  error  estimates  are 
recommended. 

statistical  error:  20  percent 
systematic  error:  negligible 
total  error:  20  percent 

These  quoted  errors  are  for  the  calibration  factor.  Cal.  To  obtain  the  error 
associated  with  measured  0^(a)  concentrations,  the  statistical  uncertainty  in 
the  measured  detector  voltage  when  viewing  the  COIL  gas  medium  must  be 
included.  Typically,  this  additional  source  of  error  is  small.  The  quoted 
error  does  not  apply  to  the  yield  of  singlet  oxygen,  as  the  uncertainty  in  any 
ground  state  oxygen  measurement  may  be  large  and  has  not  been  included. 
Finally,  yields  are  often  quoted  in  percentages  and  the  above  relative  errors 
in  tO^(a)]  should  not  be  confused  with  the  absolute  error  in  percentage 
yields. 
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6.0  DISCUSSION 

The  method  of  calibrating  O^Ca^A)  emission  detectors  by  EPR  spectroscopy  has 
been  used  by  many  research  groups  since  1978  to  support  the  development  of 
Chemical  Oxygen-Iodine  Lasers  (Refs.  4,10-12).  Typically,  these  studies 
achieved  statistical  errors  of  10  to  15  percent.  Probably  the  best  charac¬ 
terized  study  of  systematic  errors  prior  to  the  current  effort  is  documented 
in  Reference  10.  Several  discrepancies,  recommendations  for  reducing  syste¬ 
matic  errors,  and  problems  with  transfer  calibrations  are  reported  in  the 
previous  work  (Refs.  10-12).  Only  a  limited  number  of  calibrations  were 
performed  on  the  same,  unaltered  detection  system  and  thus,  little  information 
was  available  on  the  long  term  stability  and  statistical  fluctuations  in  the 
calibration  factor. 

This  study  provides  a  detailed  examination  of  the  characteristics  of  the 
calibration  apparatus  including  excited  oxygen  flow  system,  EPR  spectrometer, 
infrared  detector  system,  and  data  acquisition  system.  The  total  calibration 
error  is  largely  controlled  by  noise  in  the  observed  EPR  spectrum  of  0^(a*A) 
and  the  statistical  error  bound  is  in  agreement  with  previous  estimates. 
However,  several  improvements  in  controlling  systematic  errors  are  currently 
reported,  including:  (1)  design  of  an  optical  collection  and  detection  system 
as  a  single  unit  to  limit  the  sensitivity  of  the  calibration  to  geometrical 
and  hardware  changes,  (2)  complete  characterization  of  the  excited  oxygen  flow 
(fluid  dynamics,  kinetic  deactivation,  stability,  etc),  and  (3)  enhanced 
analysis  of  observed  EPR  spectra. 

Improvements  in  the  apparatus  and  data  analysis  have  reduced  or  characterized 
several  important  sources  of  systematic  error.  Several  advantages  are 
obtained  by  using  a  single  Infrared  detector  to  simultaneously  sample  the 
signals  from  three  viewing  volumes:  reduced  system  complexity  auid  cost, 
reduced  variability  between  detectors,  and  direct  measurement  of  relative 
signals  needed  for  an  accurate  measurement  of  kinetic  deactivation. 
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The  calibration  apparatus  has  been  designed  for  complete  characterization  of 
the  excited  oxygen  flow  including  a  direct  measurement  of  kinetic  deactivation 
between  the  diagnostic  duct  and  EPR  cavity.  The  construction  of  the  collec¬ 
tion  optics  and  the  infrared  detector  as  a  single  unit  is  absolutely  vital  in 
preventing  changes  in  detection  efficiency  due  to  hardware  configuration.  An 
analysis  of  the  EPR  signals  including  symmetry,  baseline,  truncation  far  from 
resonance,  and  line  shape  provides  increased  confidence  in  the  relative  con¬ 
centration,  r  =  [0  (a) ]/[0^(X) ] . 

No  significant  error  was  introduced  from  the  following  sources:  NO  +  0 
afterglow,  IR  detector  linearity,  EPR  linearity,  transverse  pressure 
gradients,  temporal  stability  of  0^(a)  source,  inhomogeneous  magnetic  field 
(symmetric  absorption  spectra],  and  temperature  measurements. 

Most  important  is  the  current  assessment  of  the  long  term  performance  of  the 
diagnostic.  During  a  period  of  18  months,  5  calibrations  were  performed  on 
the  same  detection  system  and  the  calibration  factors  remained  constant  to 
within  10  percent.  This  diagnostic  was  used  extensively  on  COIL  devices.  A 
complete  documentation  of  its  performance  on  the  ROTOCOIL  device  will  provided 
in  another  document.  The  consistency  of  the  present  methods  and  apparatus 
are  unparalleled  and  represent  by  far  the  most  extensive  data  base  to  evaluate 
the  performance  of  this  diagnostic. 

Several  opportunities  exist  for  increasing  the  precision  of  the  present 
diagnostic.  In  order  to  decrease  the  statistical  error,  a  higher  concentra¬ 
tion  of  singlet  oxygen  is  required  to  enhance  the  observed  EPR  spectra  of 
O^Ca).  This  cannot  be  accomplished  with  the  current  microwave  source  of 
excited  oxygen  but  might  be  achieved  with  a  higher  power,  transverse  flow 
microwave  generator  similar  to  those  currently  being  developed  to  study 
active  nitrogen.  While  a  small  chemical  generator  of  0^(a)  could  be  used, 
the  added  complexity  and  analysis  difficulties  are  considerable. 


Schlie,  L.A.,  Phillips  Laboratory,  private  communication. 
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A  calibration  of  the  magnetic  field  strength  was  not  accomplished  in  the 
present  study  and  should  be  undertaken.  Any  error  in  the  field  strength  is 
doubled  when  propagated  through  the  double  integration  and  accounts  for  the 
single  largest  potential  source  of  systematic  error. 

An  independent  determination  of  wall  deactivation  rates  would  assist  in  the 
deactivation  correction.  The  observed  dependence  of  wall  deactivation  prob¬ 
ability,  r.  on  flow  conditions  is  particularly  disturbing.  The  uncertainty  in 
kinetic  deactivation  is  small  and  only  weakly  affects  the  uncertainty  in  the 
calibration  factor.  Thus,  the  wall  deactivation  rates  are  not  particularly 
Important  to  the  calibration.  Note  that  the  pump  curve  for  this  apparatus 
provides  a  fixed  time  between  observation  points  independent  of  pressure  and 
flow  rate. 

The  NO  +  0  afterglow  was  not  a  problem  in  the  current  apparatus.  A  simple 
diagnostic  for  potential  sources  of  emission  other  O^ta^A)  is  readily  avail¬ 
able  and  is  worth  implementing  in  the  experimental  apparatus.  Specifically, 
the  simultaneous  detection  of  emission  at  1.268  /m  and  1.315  pm  is  possible 
with  two  narrow  band-pass  filters  and  a  signal  at  1.315  pm  would  indicate 
background  emission  problems. 

A  semiquantitative  assessment  of  changes  in  viewing  volume  was  made  in  the 
present  study.  A  more  complete  study,  both  analytically  and  experimentally, 
of  the  viewing  volume,  sensitivities  to  optical  focus  and  hardware  geometry, 
and  transfer  calibrations  to  different  viewing  geometries  would  be  necessary 
if  transfer  calibrations  were  necessary. 

The  prediction  of  COIL  laser  power  from  the  measurement  of  O^Ia^A)  concentra¬ 
tions  has  not  been  completely  successful.  This  may  be  due  to  several  causes: 

*  the  laser  typically  operates  near  threshold  conditions  and  small 
differences  in  the  yield  of  O^Ia^A)  produce  large  changes  in  output  power 
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•  laser  power  depends  more  strongly  on  yield  than  on  lO^Ca^A)),  and  a 
large  uncertainty  in  [0^(X)]  may  be  present 

•  a  large  uncertainty  in  several  important  energy  loss  mechanisms  such 
as  iodine  dissociation  and  nozzle  mixing  efficiency  exists 

•  COIL  devices  tend  to  operate  more  efficiently  at  higher  powers 

While  a  diagnostic  for  singlet  oxygen  is  absolutely  vital  to  characterizing 
COIL  device  performance,  its  singular  predictive  capability  is  limited. 

The  development  of  advanced  diagnostics  for  0^(X,a,b)  have  been  considered  by 
many  groups.  Two  particularly  interesting  suggestions  include  absorption  on 
allowed  ultraviolet  transitions,  and  microwave  absorption  techniques.  Absorp¬ 
tion  techniques  offer  the  distinct  advantages  of  simplicity  and  requires  no 
detailed  calibration  procedures.  However,  several  problems  exist  with  these 
approaches.  First,  accurate  determination  of  the  absorption  cross  section  is 
problematic.  Second,  pressure  broadening  effects  may  limit  the  resolution  of 
microwave  spectra.  Absorption  due  to  multiple  transitions  even  at  a  single 
wavelength  induce  difficult  deconvolution  issues.  Possibly  more  important 
than  Improved  diagnostics  for  0^(a)  is  the  development  and  implementation  of  a 
direct  0^(X)  diagnostic. 
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7.0  CONCLUSIONS 

The  precise,  absolute  detection  of  singlet  oxygen  is  a  task  requiring  great 
attention  to  detail  and  careful  experimental  methods.  Detection  by  absolute 
radiometry  calibrated  by  the  EPR  spectroscopic  method  is  currently  the  most 
reliable  method  for  quantifying  0^(a)  concentrations  for  chemical  laser 
applications.  A  calibration  apparatus  and  detector  assembly  has  been 
developed,  characterized,  and  documented  that  provides  for  reliable  imple¬ 
mentation  on  chemical  oxygen-iodine  lasers.  This  improved  apparatus  provides 
for  direct  measurement  of  kinetic  deactivation,  multiple  detections  with  a 
single  infrared  detector,  and  comprehensive  characterization  of  the 
calibration  equipment.  Concentrations  of  singlet  delta  oxygen  in  the  range 

13  17  3 

of  7  X  10  -  4  X  10  molecules/cm  have  been  measured  to  within  an  error 

bound  of  ±  12  percent.  The  precision  of  the  current  calibration  is  limited 
primarily  by  the  signal  to  noise  ratio  of  the  observed  0^(a)  EPR  absorption 
spectra.  Recommendations  for  Improving  the  precision  of  the  diagnostic  are 
presented.  Since  the  performance  of  COIL  devices  is  primarily  controlled  by 
the  yield  of  singlet  oxygen,  [O^Ca) ]/[0^(X) ] ,  the  development  of  a  diagnostic 
for  02(X),  or  yield  directly,  should  be  a  priority  for  COIL  research. 
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